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Translocation of the Global Positioning System has proved to be a highly 
accurate method of position determination for onshore and airborne navi- 
gation, but it had not been previously evaluated for nearshore hydro- 
graphic surveys. The technique of translocation for hydrographic opera- 
tions involves the simultaneous reception of signals from the GPS 
satellites by two independent receivers; one receiver onboard the survey 
vessel and one located at a known stationary position. A position 
correction is obtained at the stationary receiver, which is then applied 
to the shipboard receiver online or during postprocessing. An accuracy 
determination of translocated GPS was conducted at the Naval Postgraduate 
School, Monterey, CA, May, 1981. Two methods of positioning were used 
for comparison with GPS: 1) a least squares solution of three lines of 

position observed from three Wild T-2 theodolites, and 2) a position 
determined from ARTEMIS, a range-azimuth short range microwave position- 
ing system. Translocated GPS accuracies of 10 meters were determined. 

It is anticipated that greater accuracies will be obtained by using a 
more sophisticated receiver and more advanced processing methods. 
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I. INTRODUCTION 



Shipboard operations for the purpose of oceanographic 
research, commercial transport, tactical operations, or 
hydrographic surveying require accurate and dependable 
positioning systems. The operation of most electronic 
positioning systems is hampered by atmospheric effects on 
the propagated signal and design limitations which restrict 
operating ranges, positional accuracy and dependability. A 
universal positioning system that provides accurate 
information to the user regardless of atmospheric conditions 
and that is free of range limitations is desirable. With 
the advent of satellite positioning techniques, such a 
system is now available. 

The Navy Navigation Satellite System (NAVSAT) , or 
TRANSIT, has been operational since 1964 and has proved 
dependable and sufficiently accurate for certain users. 
Although developed for military uses, the application of 
TRANSIT in the civilian community has continued to expand at 
an exponential rate. The system, however, is limited by the 
time interval between fixes, which is 90 minutes on the 
average, and the requirement for precise knowledge of the 
ship's velocity. Although the system does not provide 
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One method of improving system accuracy is by using two 
receivers simultaneously in a technique termed 
translocation, or differential mode. The method is based on 
the assumption that signal propagation errors to two 
proximate receivers are nearly the same. One receiver is 
placed over a known position, and the variation of the 
observed position with time is determined relative to the 
absolute position. This variation is applied to the 
observed positions determined by a simultaneously operating 
receiver located at another location — on a survey vessel, 
for example. The assumption underlying this technique is 
that the distance between the two units is not sufficient to 
introduce a significant additional propagation error. The 
resulting position of the translocated receiver should be 
more accurate than the position obtained using a single 
receiver. 

For hydrographic operations, the translocated Global 
Positioning System (GPS) provides an alternative to 
presently used positioning systems. The satellites orbit at 
20,000 km, which theoretically would allow receiver 
separation of as much as 500 km without causing significant 
additional system error. With this flexibility, one 
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sufficient positional accuracy for hydrographic surveying, 

* 

it has been used for geodetic application by the Defense 
Mapping Agency (DMA) using fixed shore-based receivers to 
establish positioning networks in inaccessible areas 
worldwide. [Ref. 1] 

The Global Positioning System, or NAVSTAR, is a more 
sophisticated system designed to replace TRANSIT in the 
mid-1980s. Six satellites are currently in orbit with 
twelve more to be launched before the operational date of 
1987. The 18-satellite orbit configuration will provide a 
multi-user, passive system with 24-hour availability, global 
coverage, and possible fix intervals of less than one 
second. The operational network consists of three 

elements — the satellites, the ground-based control segment, 
operated by the Space and Missile Systems Organization 
(SAMSO) , and the user's receiver. Prototype receivers have 
been tested under optimum conditions by SAMSO at the Yuma 
Proving Ground, Arizona, on both mobile and stationary 
platforms. Results show excellent system stability and 
accuracy when compared to currently used positioning 

systems. [Ref. 2] 
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shorebased receiver could be placed at a central location 
and a survey of several hundred kilometers of coastline or 
offshore area could be completed with only one additional 
receiver on the survey vessel. [Ref. 3] The value of this 
flexibility would be realized in the significant decrease in 
survey planning, cost, and time compared to that spent 
presently in conducting hydrographic operations. 

The accuracy of the system would be sufficient for all 
survey scales. The National Ocean Survey (NOS) defines the 
required accuracy of a positioning system to be no more than 
one-third of the survey accuracy, which is 1.5 millimeters 
at the scale of the survey. [Ref. 4] Non- translocated GPS 
has been shown to be of sufficient accuracy for 1:80,000 
scale surveys and smaller. [Ref. 5] 

The object of our research is to determine whether GPS 
at the present operational level is adequate for 1:10,000 
scale surveys. To do this, a position determined by two GPS 
receivers in the translocated mode is compared to the same 
position determined by two other methods. If translocated 
GPS can be shown to achieve specified inshore survey 
accuracies, the system would provide sufficient accuracies 
at all smaller scale surveys. 
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Part of the difficulty in determining this level of 
accuracy concerns the reference system used. Host geodetic 
positions published by the National Geodetic Survey (NGS) 
are based on the North American Reference Datum of 1927 
(NAD-27) and the Clarke 1866 Ellipsoid, the best 
approximation to the shape of the earth for the continental 
United States. The satellite systems, however, use a 
reference system which is not based on a portion of the 
earth’s surface. A mass-centered ellipsoid referenced to 
the earth’s center of mass has been computed using observed 
gravity and astronomical data. The best fit to this global 
ellipsoid is the World Geodetic System of 1972 (WGS-72) , 
which is utilized by both ' TRANSIT and GPS. Transformation 
equations exist which convert NAD-27 coordinates to WGS-72 
coordinates, but uncertainty in the gravity model produces 
errors as large as 10 meters due to the transformation 
process alone. [3ef. 6] 

To establish all positions on the same reference 
ellipsoid for this research, a TRANSIT receiver was acquired 
from DHA and a first order Doppler station was established 
at the central location of the test area. A third order 
geodetic survey was conducted from the TRANSIT station to 
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establish two additional stations on the WGS-72 ellipsoid. 
Two single channel Manpack GPS receivers were obtained from 
SAHSO. Stability testing and a receiver comparison test 
were conducted prior to shipboard testing to determine a 
quantitative measure of receiver performance. One of the 
receivers was then placed over the TRANSIT station mark and 
the other aboard a survey vessel which maneuvered along 
predetermined track lines. Data were recorded via an 
interface to an HP-9825 computer and printer at each site. 
Simultaneously, two additional independent measurements were 
made of ship position — one using a range-azimuth short 
range microwave positioning system (ARTEMIS) and the other a 
position determined by a least squares solution of three 
lines of position observed from three 1-second theodolites 
set on the three pre-established WGS-72 locations. The 
comparison of the ARTEMIS and theodolite-derived ship 
positions with the translocated GPS ship position provided 
the test results. 
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II. EOOIPMENT DESCRIPTION AND OPERATIONAL CAPABILITIES 



Two positioning systems were chosen for comparison with 
the Global Positioning System: a theodolite network and a 
microwave positioning system. The theodolite network 
consisted of three shore-based units so as to provide three 
lines of position. The microwave positioning system 
provided, a range and an azimuth from a known location on 
shore. Both systems are accurate at short ranges and 
represent two contrasting methods of positioning — visual and 
electronic. The theodolite method of positioning at short 
range has been and continues to be used for hydrographic 
operations as an alternative to electronic systems when 
geometry and system accuracy limitations impose the need for 
a more flexible positioning system. The requirement for 
experienced theodolite observers as well as the 
communication and logistics difficulties involved limit the 
use of this method to special cases. 

This chapter describes the ARTEMIS short-range 
electronic positioning system, the operation of the 
theodolite network as used during the comparison test, and 
the TRANSIT satellite system. 
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ARTEMIS POSITIONING SYSTEM 



The ARTEMIS position fixing system is a short-range 
licrowave system (to 30 km), built by Christiaan 
fuygenslaboratorium B. V., in the Netherlands. It has been 
ised operationally (primarily in Europe) since 1972, when 
.he system was introduced at the International Hydrographic 
:onference in Monaco. At the present time, fifty systems 
ire in use. The system is unique in design, is highly 
iccurate at short ranges, and is particularly useful for 
small harbor surveys, river surveys and relative positioning 
lses. For this study, the ARTEMIS system, along with 
technical assistance, was obtained from the American 
listributor, MARINAV Corporation of Houston, TEXAS. 

ARTEMIS is a range-azimuth positioning system which 
employs a unique tracking method. A stationary or fixed 
unit is placed over a known position on shore. Using a 
telescope that is mounted on the antenna, the antenna is 
sighted onto another known point and the known azimuth 
between the two locations is entered into the system via a 
digital display on the fixed unit. The unit is then locked 
into place for operation. A mobile unit is placed on the 
vessel whose position is desired. 
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The mobile and fixed antennas resemble radar antennas 



which, while operating, track each other so that the two 
antennas are always parallel to each other and perpendicular 
to the line of position between them. A microwave link is 
established between the two antennas during periods of 
measurement and information is relayed between the mobile 
and fixed units. (A voice link is also available during 
operation.) The fixed unit provides the direction and the 
azimuth of the mobile unit with respect to the point sighted 
on by the telescope. At the same time, the mobile unit 
transmits a coded interruption of the signal to the fixed 
unit, which replies with the same interruption. The time 
difference between the transmission and reception at the 
mobile unit determines the distance between the two 
a ntennas. 

Distance accuracy is determined by the number of 
measurements averaged over time and is available in two 
modes — a dynamic mode which averages 1000 measurements or a 
static mode which measures 10,000 measurements per displayed 
distance. The azimuth and distance information is displayed 
by the mobile unit so that personnel are not required at the 
stationary unit once operation has begun. The system can be 
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operated in long or short-range modes, which determine the 
amount of power output for the system. 

The microwave frequencies used are 9.2 and 9.3 GHz. 
Problems arising with signal interference, observed with all 
microwave positioning systems, therefore can occur with 
ARTEMIS. The regions of range holes, or multipath 

interference, can be computed for up to 10 Jem range by prior 
knowledge of the antenna elevations above sea level and the 
operating distance from the shore-based antenna. [Ref. 7] 

Each antenna is designed with a 22° vertical beam width 
and a 2° horizontal beam width. A power source of 22 volts 
DC per unit is required, with- an average current load of 2.5 
amps. 

The accuracy of the system depends on both the distance 
and the angle measured. The angle measurement has a 2 a 
probability of measuring within 2* of arc (0.033° ) . The 
distance measurement has a 2 ° probability that the mean 
distance value will vary within + 1.5 meters. 

One of the advantages of ARTS MIS over other electronic 
systems is that the two lines of position defining the fix 
always intersect at 90° . A unique advantage is that the 



coverage around the fixed unit can be circular. 



with a 



maximum range cf operation to 30 km, dependent on antenna 
heights. The short range applications of the ARTEMIS fill a 
need not provided by many currently used positioning 
systems. 

For this study, the ARTEMIS was used in conjunction with 
the positions derived from the theodolite network to 
determine the best estimate of the ship's position. The 
accuracy of the ARTEMIS was more than adequate for the 
ranges encountered during the test, and the ease of use and 
operational versatility of the system made its use very 
advantageous. 

B. THEODOLITE NETWORK 

A network of three Wild T-2 theodolites was used as the 
control positioning system. The theodolite is widely used 
for surveying as a means of measuring angles between two 
observational locations — either in a horizontal plane for 
traverse or triangulation surveys, or a vertical plane for 
astronomical observations. 

The use of three theodolites as a short-range 
positioning system was an accepted method of positioning 
prior to development of electronic distance measuring 
systems. There remain disadvantages when using electronic 
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positioning systems for hydrographic operations at short 
ranges (1-2 km). Most short-range electronic systems are 
designed for measurement of position at line-of-sight 
ranges. This region includes a large portion of inshore 
hydrographic surveys. Inshore positioning involves more 
thorough planning to obtain the required accuracy for the 
scale of the survey. For example, in a harbor or cove it 
may be difficult to place two range-range electronic units 
so that the area to be surveyed is not within the region of 
positioning system uncertainty, or 30° (radially) on each 
side of the baseline between the two units. To reposition 
the electronic units as often as is required to obtain 
optimum geometry is time-consuming and requires manpower and 
extensive planning, and can result in operational delays. 
Interference of the signal due to reflection from a calm sea 
surface, or multipat hing, is also a problem at short ranges 
for microwave systems. 

Theodolites can be set at any intervisible location with 
preplanned geometry and can be used for the entire region to 
be surveyed. This method requires less expensive equipment 
(three theodolites) , more manpower (a recorder and an 
observer at each site) , communication between the three 
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locations (via portable radio, for example) and reasonable 
weather conditions. The basis of this procedure is that 
three angles observed from known locations to a 
predetermined point on the vessel constitute three lines of 
position. After removing all possible systematic errors and 
through a least squares solution of the three lines of 
position, the position of the vessel is determined. A major 
disadvantage is that observational errors are not determined 
in real time because the data cannot be processed 
immediately. 

The basic procedure for the method requires 
predetermined sites which are known relative to each other. 
This may be achieved by using published geographic positions 
or by conducting a local plane survey to establish new 
positions using preestabiished stations nearby. Either way, 
all three locations must be based on the same datum and 
preferably permanently positioned by a disk or a station 
mark. 

Each theodolite was aligned using a plumb bob or optical 
plummet over one of the known station marks or disks. By a 
predetermined scheme, each site was marked so as to be 
clearly visible to the other sites (lights were used for 
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this study) . The telescope of each theodolite was sighted 
onto one of the other sites and the horizontal circle set to 
read an arbitrary value close to 0°. The crosshairs in the 
telescope were set to intersect the mark. When the 
theodolite was exactly sighted, the horizontal circle was 
read and recorded. (At this time the lights were turned 
off.) The operation began when all three sites completed 
this step, called the "initialling procedure", or "initial 
pointing". The instrument was sighted on a point on the 
survey vessel that was visible to all operators. The ship 
was tracked by the operators by using the micrometer wheels 
on the instrument until a "mark" was called over the radio 
by the recorder at the control theodolite location. At that 
instant, the operator stopped tracking and the horizontal 
circle was read and recorded. This observation constituted 
a fix. After a series of fixes were recorded, or when the 
survey vessel reached a predetermined location, a set of 
fixes was completed. The lights were then reactivated at 
each site, each observer intersected the original mark, the 
horizontal circle was again read and the "final pointing" 
was recorded. If necessary, the instrument was relevelled. 
The process was then repeated by starting with the 
initialling procedure. 
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It was assumed that the theodolites were in good 
operating condition. Collimation error had previously been 
determined. The observers were kept at the same locations 
throughout the test to minimize the error due to individual 
observer variations. Since only one portion of the 
horizontal circle was used, and no reverse readings 
observed, a systematic error was introduced. This was 
assumed to be minimized during computation of the vessel 
position. 

The observations were corrected for collimation error 
and for pointing error by evenly distributing the difference 
between the initial and final pointing values. The 
resulting angles were processed using an intersection 
computation with a least squares fit applied to optimize the 
result. The computations for this study were processed 
using an HP-9815 computer and software prepared by Cdr 
Ludvig Pfeifer, NOA A. The resulting position was used as 
the best estimate of the position of the ship. 

C. TRANSIT POSITIONING SYST2H 

The TRANSIT, or Navy Navigation Satellite System, is the 
most widely utilized satellite positioning system at the 
present time. The network of satellites and ground tracking 
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stations has been in operation since 1964 and is becoming 
widely used by both the military and civilian communities. 
It is expected to continue operation until at least 1995, 
when the Global Positioning System is expected to replace 
it. 

The system consists of a series of T3ANSIT satellites 
(fig. 2.1.) in fixed polar orbits of 1095 km elevation, 
forming a network within which the earth rotates (fig. 2. 2.). 
The orbits are spaced so that satellite passes occur, for a 
given ground location, at intervals of between 35 and 100 
minutes. Each satellite travels an observable distance of 
between 4400 and 7000 km per pass, providing a sufficient 
time interval and baseline over which the user can record 
data. 

The idea of positioning an observer on the earth's 
surface by using the satellite as a reference system 
originated from observation of Sputnik 1, launched by the 
Soviet Onion in 1957. By recording the change in frequency 

of the satellite signal as the satellite passed overhead, a 
Doppler shift could be measured. The Doppler shift gives a 
spatial relationship as a function of time between the 
observer and the satellite. If the position of the 
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Fig. 2.1. TRANSIT Satellite 
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satellite is known at the time of observation, the 
observer’s position can be determined. 

In the TRANSIT system, the satellite positions are 
precisely determined by tracking stations located at 
Prospect Harbor, Maine; Rosemount, Minnesota; and flahiawa, 
Hawaii. The tracking information (in the form of the 
Doppler frequency shift as a function of time) is relayed to 
the computing center at the Naval Astronautics Group 
Headquarters, Pt. Mugu, California. Here the actual 
satellite orbits are computed. This information is used to 
predict orbit information. Parameters of this predicted 
orbit are relayed to two injection stations at Rosemount and 
Pt. Mugu. During the next satellite pass, the information is 
relayed to the satellite as part of the navigation message, 
with new updates occurring every 12 hours. 

The satellites therefore contain predicted orbit 
information which is received by the user via a satellite 
receiver. Several receivers are available with complexity 
and positional accuracy dependent upon the needs of the 
user. A computer is normally required to determine the 
user's mobile or fixed position in real time from a single 
satellite pass. 
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For many users having a dual channel receiver and known 
velocity, this method of positioning provides sufficient 
accuracy with typical error of 27 to 37 meters root sum 
square. Optimum accuracy is obtained, however, if the 
receiver position is fixed, and several consecutive passes 
of the same satellite are recorded. During each pass, 
lasting from 10 to 16 minutes as the satellite travels from 
one horizon to the other, many Doppler counts (observations 
of shift in frequency as a function of time) are recorded. 
After a predetermined number of passes are stored, the data 
can be forwarded to the Defense Mapping Agency in Brookmont, 
Maryland. There the precise ephemeris data from the ground 
tracking system is applied to the data from the recorded 
passes to provide the optimum solution. The resulting 
position has a typical error of 6.3 meters root mean square 
for a series of passes, with a 3-dimensional result of 1.5 
meters per axis repeatability after 25 precise ephemeris 
passes. This method is used to determine positions in 
remote areas by DMA geodetic survey groups using the 
Magnavox AN/P5R-14 Geoceiver (fig. 2. 3.). [Ref. 8] 

The TRANSIT system uses the World Geodetic System of 
1972 (WGS-72) reference ellipsoid which is based on the 
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Universal Space Rectangular Coordinate System — an earth, 
mass-centered coordinate system. The WGS-72 ellipsoid 
describes the best fit of the entire earth to a mathematical 
model. The satellite orbits are computed using this 

geopotential model, the accuracy of which is based on 
gravitational and astronomical data. 

For the present study, a Magnavox AN/PRR-14 Geoceiver 
was acquired from the Department of Satellite Geophysics, 
Satellite Tracking Branch, DMA. The receiver was placed 
over a fixed point and 40 passes of one satellite were 
observed. The procedures described previously for an 
optimum solution of point positioning were followed. The 
resulting latitude, longitude, and height relative to the 
MGS-72 ellipsoid are coordinates of a first order TRANSIT 
position as defined by DMA standards. This position was 
used as the basis for horizontal control for the study (see 
test procedures) . 
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GPS SYSTEM DE SCSI PTION AND OPERATIONAL CAPABILITIES 



A. SYSTEM DESCRIPTION 

NAVSTAR , or the Global Positioning System, was designed 
for maximum operational capability — high accuracy, 
twenty-four hour access, passive user, and world-wide 
coverage. A brief description of the GPS, and the 
satellites (NAVSTARs), will serve as a background for the 
more in-depth discussion immediately following. 

GPS navigation is based on the measurement of four 
ranges, or pseudo-ranges, from each of four satellites with 
known, or predicted, positions. Each range determination is 
computed using the propagation velocity of the signal 
multiplied by the sa tellite-to-user travel time. Satellite 
transmission time is included as part of the navigation 
signal message. From the four known ranges, the user 
position and user clock error are computed on a three-axis, 
orthogonal coordinate system. This computation is then used 
to transform the user position onto the WGS-72 ellipsoid. 

The remainder of this chapter will provide the 
background for GPS development and the three major system 
segments of the GPS: the space system segment, the control 
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segment, and the user system segment. 



Included in the user 



system segment will be a discussion of the Manpack receiver 
as developed by Texas Instruments Corporation used for this 
study. 

With the successful operation of the TRANSIT system in 
the early 1960's, the Department of Defense (DOD) began to 
consider the specifications of, and the technology required 
for, the development and implementation of a second 
generation system. On 17 April, 1973, the Deputy Secretary 
of Defense stated that a Defense Navigation Satellite System 
was to be developed. [Ref. 9] Using the new synthetic 
oscillator technology, and the Department of Defense surface 
and nearsurface navigation objectives and requirements, a 
DOD interdepartmental task force developed the GPS as the 
second generation satellite navigation system. With the 
United States Air Force (USAF) as the managing agency, the 
GPS became a joint service program, with representatives 
from the DMA, Army, Navy, Marines, Coast Guard, and the 
North Atlantic Treaty Organization. [Ref. 10] The Joint 
Programs Office was established within the Air Force Space 
and Missile Systems Organization to design, develop, and 
implement the GPS program. 
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GPS was conceived and designed to provide continuous, 
real-time, world-wide navigation coverage for subsurface, 
surface and aearsurface operational vessels, vehicles, and 
aircraft. Accuracy of better than 10 meters was desired 
under all conditions, with navigational capabilities to 
update and interface with other forms of existing and 
projected navigational systems. Another requirement of the 
GPS satellites was nearly trouble-free and simple user 
operation, even in low signal-to-noise environments. An 
encoded signal structure was required for the satellite 
system for security reasons. In addition, the signal 
transmitted by the HAVSTAR satellites had to include other 
operational information in conjunction with the navigational 
message. To provide for navigational requirements in event 
of hostilities, the GPS was designed to permit the use of 
redundant satellites. 

To achieve these objectives and requirements, the GPS 
program was organized according to three developmental 
tasks: concept validation, full-scale engineering 
development and systems tests, and production and operation. 
Starting in 1974, the concept validation task was initiated 
with the merger of the TISATION series satellite program and 
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the US AF System 21B (a highly accurate three-dimensional 
navigation system) program. Tests with the Navigation 
Technology Satellite, NTS-1, utilizing two rubidium clocks 
(accuracy 5-10 parts per 10 13 ) and an orbit 7500 nautical 
miles (nm) above the surface and three ground based 
receivers (called inverted range) provided concept 
verification. [Ref. 11] Additional tests were conducted 
with an improved NTS-2 satellite positioned at an altitude 
of 10,980 nm and containing two cesium clocks (accuracies on 
the order of 1 part per 10 13 ). [Ref. 12] Orbital height of 
10,980 nm permitted the evaluation of the effect of the 
planned GPS constellation altitude on navigational 
accuracies using the improved cesium frequency standards. 
In 1980, NTS -3 was launched and the Advanced Developmental 
Model (ADM) hydrogen maser frequency standard evaluated for 
use on the GPS (providing an expected frequency accuracy of 
1 part per 10 13 ). [Ref. 13] The results obtained from the 
NTS-3 satellite have not yet been released. Since 1978, 
additional satellites have been orbited to provide the 
six-satellite constellation required to fully evaluate the 
first development task — concept validation. Phase I. 
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During Phase I, user requirements were consolidated to 
allow the development and construction of four basic types 
of test receivers: continuous tracking receivers, 
sequential tracking receivers, lew cost receivers, and the 
Hanpack receivers. [Ref. 14] These receiver types were 
then constructed for Phase I testing by various contractors 
and their performances evaluated aboard helicopters, 
aircraft, ships, vehicles and personnel at the Joint 
Programs Office's testing facilities. Phase I, concept 
validation, was completed in late 1979 and Phase II testing 
was begun in early 1980 with full scale engineering, 
developmental, and systems tests. 

Phase II tasks, finishing the operational and control 
segments of the testing, are scheduled for completion in 
mid-1983. The operational control station will be relocated 
at a central continental United States site during this 
phase. Additionally, all satellite monitor receivers at the 
remote monitoring stations will be updated to operational 
status. User equipment and receiver selection will be made 
in two increments. Four contractors will be asked to 
compete for the second increment of the receiver selection 
process by providing upgraded and refined user receivers for 



42 



extensive testing 


and 


eval uation. 


Two 


of the 


four 


contractors will 


then 


be selected 


to 


continue 


the 



competition and will provide additional receiver design, 
refinement and modification, leading to prototype production 
and extensive field testing. 

When Phase II is completed in mid-1983. Phase III, 
production, will commence with the selection of one of the 
two contractors to produce the user receivers. Manufacture 
of these receivers will commence shortly after selection, 
and commercial availability of receivers is expected in 
1985. [Ref. 15] The GPS should be fully operational with 
world-wide coverage using a constellation containing a total 
of eighteen satellites in 1988. [Ref. 16, 173 

The GPS is considered to be composed of three 
inter-related system segments: the space system segment, 
the control segment, and the user system segment. Of these, 
the space system segment may be considered the most 
important since the other two require the space segment to 
be operational at all times. Operationally, the space 
system segment consists of 18 satellites placed in orbits at 
altitudes of 10,898 nm and having orbital periods of about 
12 hours. [Ref. 18] 
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The original plan for the satellite constellation 
consisted of 24 satellites in three planes of 8 satellites 
each. These three planes were to be placed 45° of longitude 
apart and at an inclination of 63 to the equator. 
[Ref. 19] As a result of 1980 budget cutbacks the space 
system segment has been reduced to a total of 18 satellites. 
[Ref. 20] Several different types of 18 satellite 
constellations have been proposed, but all have drawbacks. 
One which appears the most promising is the 'Honuniform 18* 
constellation containing 6 planes, each of which contains 3 
satellites spaced 120° apart (fig. 3.1.). The areas of user 
outages (fig. 3. 2.) are a disadvantage of this configuration. 
The user outages are indicated by the darker areas shown on 
Figure 3.2. These outage areas will appear to rotate about 
the globe at the indicated latitudes, and positioning with 
the GPS will be intermittent, or of decreased accuracy for 
users at these latitudes. The final constellation 
configuration is to be selected early in 1982. [Ref. 2 VJ 

Sach NAVSTAR satellite has the following 
characteristics: weight of 982 pounds, length of 17.5 feet 
from tip to tip of the 5 square solar panels, three atomic 
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frequency standards (clocks) , and three nickel-cadmium 

i 

batteries for operations in darkness and during peak load 
periods. Stabilization of each NAVSTAR is maintained by 
four skewed reaction wheels which produce three-axis 
stabilization. NAVSTAR station-keeping and momentum damping 
capability is provided by a hydrazine propulsion system 
within each satellite. Life expectancy of each unit is five 
years, and expendables are expected to last for up to 7 
years with normal usage. 

Figure 3.3. illustrates the parts of each NAVSTAR. A 
twelve element shaped-beam helix antenna transmits the 
navigation message on two L band frequencies to user system 
segments. Simultaneously it receives satellite status, 
clock corrections, ionospheric data, and ephemeris constants 
from the control segment on two S band frequencies. 
(Ref. 22] The two coherent L band frequencies are centered 
on 1227.6 MHz (LI), and 1575.4 MHz (L2) . [Ref. 23] 
(Ref. 24] Both frequencies are necessary to determine the 
total ionospheric delay correction for signal propagation 
time (Table I ) . Brrors produced by the ionospheric delay 
in signal velocity contribute most of the error in the 
pseudo-range. By utilizing the LI and L2 frequencies and 
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Fig. 3.3. NAVSTAR Satellite 
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TABLE I. Ionospheric Group Delay as a Function of Frequency 



'C. = 

GD 






A 

f 



2 



(neglecting third order and higher terms which contribute minimal 
error) 



where A = ionospheric delay constant of conditions 
R = true range 
C = speed of light 
f = carrier frequency 

B « (averaged earth magnetic field strength along 
the path) 

^Gq = 9 r oup delay between points G and D 



Term Neglected: 
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C 



Range Error at 
f - 1 .5 GHz 

i 1 inch 
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the following algorithm, the error is decreased but not 
totally eliminated. [Ref. 25] 



T ' - T = A T 

GDLi GDL 2 




where A = ionospheric constant and T , T are the 
respective ionospheric group delays at frequencies LI and 
L2. Both frequencies are also modulated to transmit the GPS 
time codes and the navigation message. 

GPS time is transmitted using two pseudo-random noise 
chip codes — the Coarse Acquisition (C/A) and the Precision 
(P) codes. Modulation on the Li frequency carries both the 
C/A and P codes in phase quadrature, while the L2 frequency 
carries only the P code. These codes identify each acquired 
satellite by matching the unique pseudo-random noise code 
pattern generated by each satellite with similar user 
generated codes. 
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Both codes are required to measure the navigation signal 
propagation time. The phase shift is measured and used to 
match the user generated signal with the incoming satellite 
signal. [Bef. 26] The C/A code has a repetition period of 
one millisecond and has a short code stream, transmitting 
1.023 million bits per second (bps). This code is easily 
acquired and matched with the user generated pseudo-random 
noise code but only provides a coarse time signal for 
computing the pseudo-range. However, the P code has a long 
code stream, repeating only every seven days, and it is 
transmitted at the higher rate of 10.23 million bps. With 
this higher data transmission rate and long repetition 
period, a more accurate time signal is available for 
pseudo-range determination. However, the long period and 
high data rate makes the P code difficult to match with the 
user-generated pseudo- random noise code pattern unless the 
receiver has a highly accurate time standard, and the 
approximate receiver location is known. In lieu of this, 
the usual method of acquiring the P code is by matching the 
C/A code of the desired NAVSTAR, then using the Handover 
Word (HOW) , to locate the correct P code sequence within the 
7-day pseudo-random noise code. GPS time is transmitted 
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with the HOW every six seconds and provides very accurate 
GPS time to the user. 3oth of these codes provide the 
proper time reference for computation of the pseudo-range. 
The C/A code provides a less accurate time reference than 
the P code. The selection of the code depends only on the 
navigational accuracy desired. 

The navigation message is also modulated on frequencies 
LI and L2 but at a much slower rate (only 50 bps) than 
either the C/A or P code. On both frequencies, the 
navigation message is identical and consists of 5 subsets of 
6 seconds duration each. Each navigation message begins 
with the telemetry message (TLM) and the HOW message, which 
are used to transfer from the C/A to ? codes. The remainder 
of the entire navigation message allocated to each 
subsection is transmitted after the TLH and HOW messages. 
Information in the navigation message includes approximate 
satellite ephemerides, status of other NA7STA3 satellites, 
parameters for clock corrections, atmospheric signal delay 
corrections, satellite performance status, momentum dump 
status since last upload, time since last upload, almanac 
and identification codes for all satellites, and provisions 
for inclusion of any special or important messages from the 
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control system segment to the user. [Ref. 27, 283 
Navigation is thus achievable by combining the information 
transmitted by the space system segment in the form of the 
C/A and P codes and the navigation message. 

Accuracy and proper function of GPS are governed by the 
second segment — the control segment. The control segment 
monitors and uploads the individual NAVSTAS satellites, 
checks operational characteristics, determines the ephemeris 
and almanac data and corrections required, computes the time 
delay and clock corrections for each satellite atomic 
standard, determines the atmospheric delay corrections for 
pseudo-range computations and includes special messages in 
the body of the navigational message. These functions are 
under the command and operational control of the Master 
Control Station, currently located at Vandenberg AF3 , 
California. Four unmanned remote monitor stations, located 
in Hawaii, Guam, Elmendorf AFB, Alaska, and Vandenberg AFB , 
collect the pseudo-ranges to each satellite, the change in 
pseudo-range of each satellite signal, local meteorological 
data, and the remote monitor station atomic standard 
parameters. Each set of data is transmitted, on request, to 
the Master Control Station, where all the data is combined 
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and the satellite upload corrections and parameters are 
computed. These corrections and parameters are then relayed 
to the Opload Station, Vandenberg AFB, for transmittal to 
the respective NAVSTAR satellites. The uploads occur when 
the satellites are first visible to the Upload Station and 
when the Oser Range Error exceeds 4 meters, as computed by 
the Master Control station. [Ref. 29, 30, 31H 

The las-t segment, the user system segment, combines the 
information provided to the space system segment (by the 
control segment) with the information received from the 
space system segment to obtain the pseudo-ranges from four 
MAVSTAR satellites. These ranges are computed using the 
following algorithms: 

B 'i = ' cAt A. 

where R. = true slant range to satellite i, c = speed of 
light in a vacuum, tp = time of received GPS signal by user 
receiver (assumed to be simultaneous from all satellites) , 
t^ = time of transmission at satellite i, and 
At^ = atmospheric propagation delay time from satellite i 
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(fig. 3. 4.)* Assuming a user clock error, the true range can 
be expressed in terms of a pseudo-range: 

R p = R + cAt A ♦ C(Atp- A t s ) 
l i ' l 

where Rp. = pseudo-range including clock error, At^ = time 
delay in the user receiver cf the received signal, and 
At$. = time delay in satellite transmission from satellite 

i. 

To determine the slant ranges: 

R • = n|(x s -x) 2 + (y s -y) 2 ♦ (z. -z) 2 

i i 5 i 

R i =J(X S _-X) 2 + (y s . - y) 2 +(z s -z) 2 +cAt A> +c( Atp -At $ . ) 

The following quantities are known: x , y ,2 obtained 

from navigation ephemeris data; A t A computed using the 

l 

ionospheric delay correction obtained from observing 
frequencies LI and L2 (tropospheric delay neglected) ; its. 
obtained from the navigation message as the satellite clock 
error. Thus, there are four unknowns: x, y, z (user 
position coordinates) and Atp ( user clock error or bias). 
To solve this set of equations, four satellites must be 
observed simultaneously. [Ref. 32] 
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Should only three satellites be observable, the altitude 
may be assumed and the solutions can be computed for the 
other two location variables and the user clock error. 
Solutions of this type (altitude hold) significantly degrade 
position accuracy. 

The above description is over-simplified since the 
relativistic effects of NAVSTAR satellite velocity, as 
observed by the user, are neglected in terms of the 
satellite time delay correction and the atmospheric 
propagation delay correction. The actual user velocity, in 
terms of velocity in the x, y, and z directions, is 
calculated using the Doppler shifts of the two frequencies 
Li and L2 and the estimated satellite position taken from 
the ephemeris data. With the relativistic effects as an 
integral part of the data set, solution of the navigation 
position problem becomes more complex. Phase I receiver 
manufacturers have included them in the solution of the 
problem by several different methods. Methods of solution 
depend on the usage of the receiver (whether static, high or 
low dynamic) and the expected cost/benefit ratio for the 
receiver (in terms of the total cost, weight and size) . 
Current state-of-the-art design technology and navigational 
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satellite vehicles 




Fig. 3.4. Pseudo-Range Description 



57 



algorithms, developed since Phase I, s#ill be included in the 
Phase II receivers and will significantly improve accuracy 
of future receivers. 

B. MANPACK RECEIVER 

The Phase I receiver used for this research was the 
Texas Instruments GPS Manpack, also known as the Manpack 
Vehicular User Equipment (MVUE) . MVUE is a single channel 
low dynamic receiver, usable for receiver velocities from 0 
to 25 meters per second. [Bef. 33] System configuration is 
illustrated by figure 3.5. Composing the system are the 
aanpack control data unit (CDU) , antenna and vehicle mount, 
and vehicle power converter. The Manpack is the receiver 
unit (fig. 3.6.) with provisions to connect an external data 
output, through the MVUE Instrumentation System (MIS) 
interface, figure 3.7. Facilities for an internal 
preamplified antenna mount, or an external preamplified 
antenna, and an attachment for either an external power 
source or the integral battery pack are also shown (figure 
3.7). Together, the MVUE and the SBP/9900 microprocessor 
data processing unit have storage capacity for 47,360 
sixteen bit words. This storage capacity allows 
interrogation of NAVSTAR satellites sequentially, with the 
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Fig. 3.5. Global Positioning System MVUE Receiver 
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receiver computing a position fix every 24 seconds. The 
relativistic velocity effects on the received signals are 
determined by matching the incoming frequencies with the 
expected frequencies. The user velocity relative to the 
satellite is determined from the observed frequency shifts. 
The nominal frequency shift for each satellite transmission 
is computed using the simultaneously transmitted ephemeris 
data. Total weight of this receiver unit is 33.5 pounds, 
and it has a volume of 17,958 cubic centimeters. 

The second element of the HVUE, the CDU (fig. 3.8.), is 
the input/output device for the non-automated Manpack user. 
Inputs are entered by keying the pressure pads set in six 
rows of four pads per row. To obtain a navigational 
position, the approximate user latitude, longitude, and 
Greenwich Mean Time (GMT) are keyed. After the warmup 
period of 3 to 5 minutes, the satellite signals are 
received, and the CDO will display the actual position. 
Additional input may include NAV STAS satellite selection, 
user altitude, and other position datums. Output options 
include position, GMT time, and system error messages. 
Other inputs and outputs are listed in the MVUE manual. 
(Ref. 34] 
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Fig. 3.3. MVUE CDU Keyboard/Display Layout 
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Automated users may elect output either to a radio 
interface or through the MIS interface. For this research, 
a special interface unit, designed and constructed at the 
Office of Testing and Evaluation of Systems, NOAA , was used 
to link the MIS interface board within the receiver with a 
Hewlett Packard 9825T computer. With this interface, user 
latitude, longitude and GPS time were logged onto magnetic 
tape for future processing. Figure 3.9. and Table II are 
the engineering documents used in construction of the MVtJE 
interface unit. 

Other MVOE elements are the antenna and the antenna 
external preamplifier vehicle mount. Located within the 
antenna are two antenna subassemblies associated with 
frequencies LI and L2. Both of the subassemblies are 
connected to the external preamplifier with two lengths of 
coaxial cable, enclosed in a flexible conduit. The entire 
unit, antenna and preamplifier, is connected to the Manpack 
receiver with a power cable and a coaxial antenna cable. 
Cable length is determined by the location of the external 
antenna but is limited by the minimum strength required by 
the Manpack receiver. The combined antenna and preamplifier 
unit is 36 inches long, and it weighs 1.9 pounds. 
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Circuit Diagram of Receiver Interface 



TABLE II. MANPACK Interface Wiring Code 
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The power converter must supply 24 volts DC to the 
Manpack receiver. During this research, a power supply 
which converted 120 volts AC to 24 volts DC was used. This 
power supply was considered more stable than the power 
converter provided with the MVOE, and the voltage output and 
power were controlled and monitored during the testing 
period. 

C. DIFFERENTIAL GPS 

As mentioned earlier in this chapter, atmospheric 
propagation delays are the major source of error in the 
navigational equations because ionospheric propagation is 
poorly understood. The tropospheric propagation delays are 
neglected in the navigation equations, and the resultant 
error is dependent on local temperature, pressure, humidity, 
and NAVSTAR elevation and azimuth. Tropospheric corrections 
are usually neglected if the observed satellite is more than 
5° above the local horizon. [Ref. 35] 

Additional error sources are receiver and satellite 
clock perturbations, lengths of signal paths from different 
satellites, multipathing delays, satellite ephemeris errors, 
and receiver noise. Combined, these factors lead to errors 
of 3.6 to 6.3 meters, as table III illustrates. [Ref. 36] 
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TABLE 1 1 1. RANGE ERROR BUDGET 



Uncorrected Error Source 



User Equivalent Ranae Error, 

TT^T ' 

feet meters 



SV Clock Errors * 

Ephemeris Errors 5.0 

Atmospheric Delays 8.0 - 17.0 

Group Delay (SV Equipment) 3.3 

Multipath 4.0 - 9.0 



1 .5 

2.4 - 5.2 

1 .0 



1.2 - 2.7 



Receiver Noise and Resolution 
Vehicle Dynamics 



Root Sum Square 11.8 - 20.7 3.6 - 5.3 

*Two hours after update 
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One method used to reduce the magnitude of this error is 
that cf differential, or translocated, positioning. To use 
this method, one receiver is positioned over a known 

geographic position. Pseudo-ranges to each satellite are 
determined, and a differential correction is obtained for 
each pseudo-range. These corrections are then applied to 
the satellite pseudo-range observed at a second mobile 
receiver. aith the differential corrections applied, the 
net error of the mobile GPS receiver is significantly 
reduced. [Ref. 37] 
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IV. TEST PROCEDURES 



A. POSITIONING NETWORK 

One of the considerations in planning this test was to 
examine the various methods of positioning the shore-based 
equipment. It was desirable to use closely spaced (1-2 km) 
locations within the immediate surroundings of Monterey 
harbor for reasons of ready access by personnel during 
installation of equipment and for optimum accuracy of the 
positioning systems chosen for control. As mentioned 
previously, the positioning systems chosen for the 
comparison (the three theodolites and the ARTEMIS) give best 
results at relatively short ranges. 

Two predetermined station marks were selected for the 
test — 1) USE MON 1978, a bronze disk set by the Army Corps 
of Engineers, located on a sand dune on the Naval 
Postgraduate School property 400 meters from the ocean at an 
elevation of 40 meters above sea level and 2) MUSSEL 1932, a 
bronze disk set by the Coast and Geodetic Survey in a rock 
outcrop 7 meters above sea level at the shoreline on the 
property of the Hopkins Marine Laboratory of Stanford 
University. The third station mark was set in the center of 
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the Monterey Coast Guard Pier for purposes of the test. 
This was a DMA bronze disk stamped GE0CEI7ER STA 31370. 
This location was chosen foe access to power facilities 
available at the pier and because of central location to the 
test area. 

The position of the Geoceiver station was established 
using a Magnavox Geoceiver (Model AN/PRR 14) obtained from 
the Department of Satellite Geophysics, Satellite Tracking 
Branch, DMA/Hydrographic Topographic Center. The equipment 
was operated in accordance with DMA requirements to obtain a 
documented first order station. The data were submitted to, 
and processed by, DMA/HTC in Brookmont, Maryland, using 
precise ephemeris information to provide the WGS-72 
position. Michael Ellett of DMA provided the technical 
expertise for acquisition of the data and maintenance of the 
geoceiver during the test period. 

A third order. Class I, closed traverse was conducted 
following the standards of the National Geodetic Survey 
(NGS) , using a Wild T-2 theodolite to measure angles and a 
Motorola MRA V Tellurometer to measure distances. The 
TRANSIT position computed by DMA/HTC was used for the 
initial position of the survey, and a Polaris observation 
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using the theodolite provided the initial azimuth (fig. 
4.1.). A total closure error in distance of 0.061 meters 
and angle of 5.4" was achieved. The positions determined 
from the traverse and used for the control network are as 
follows: 

USE HON 36° 36 ' 08 .424"N 1 2 1° 52 • 4 0 . 30 1" W 

MUSSEL 36° 37' 21 .882"N 1 2 1° 54 ' 16 . 048"W 

GEOCEIVER STA 31370 36° 36' 36 .246"N 1 21° 53' 29 .693"W 

B. POSITIONING EQUIPMENT 

As previously described, three positioning systems were 
used during the test — the three theodolites, the GPS, and 
the ARTEMIS. In addition, the geoceiver was set to acquire 
data during the shipboard portion of the test. During the 
time of day that the GPS satellites were available (Table 
IV) only one or two observations per day of the same 
TRANSIT satellite were possible. The geoceiver data were 
recorded but were not used for this study. 
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TABLE IV. SATELLITE OBSERVATION TIMES 
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1 . Theodolites 



This section describes the placement and operation 
of the equipment during the test. One theodolite was set on 
a tripod at each of the three stations, and each was 
occupied by an observer and a recorder. Communication 
between stations was maintained via portable radios. 

Because tests could be conducted only between 10 PM 
and 2 AM local mean time (LMT) due to GPS availability 
times, it was necessary to plan a systems of lights for the 
stations. The three theodolite sites were made visible by 
placing battery powered lights directly over the station 
marks. In addition, battery operated light packs for the 
theodolites were used. All shorebased lights and light 
packs for the test were provided by the Pacific Marine 
Center, NOAA. 

The ship system was designed with an omnidirectional 
pair of lights anchored around the support pipe for the GPS 
antenna. Also, the equipment shelter on the boat deck of 
the ship was lighted by spotlights so that the GPS antenna 
was easily seen by the shorebased theodolite operators. 

During the test, each theodolite observer initialled 
on one of the other stations and recorded the reading on a 
log sheet. The theodolite operators then intersected the 
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lights on the ship GPS antenna using the theodolite 
telescope and tracked the ship until a 'nark* was called 
over the radio from the personnel aboard the ship. The 
angle measured at this time was recorded on the log sheet as 
a fix. The time interval between fixes was approximately 
one minute. Each series of approximately twenty fixes 
constituted a set. At the end of each set, a final reading 
was made of the angle to the reference station. The fixes 
were numbered consecutively throughout the test, and the 
Greenwich Mean Time (GMT) of each fix was recorded on the 
ship log sheets. The length of the sets depended on the 
ship position and heading. 

To correct the observed theodolite angles, the mean 
of the final pointing and the initial pointing was 
determined for each set and for each station. This value 
was assumed to be a constant throughout the set and was 
removed from each angle. A collimation error for each 
theodolite was also removed from each theodolite 
observation. The same observers and theodolites were used 
nightly at all stations except MUSSEL, where two different 
individuals observed over the test week. This continuity of 
personnel and equipment was designed to reduce operator and 
instrument error to a minimum. 
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2. GPS and ARTEMIS 



To colocate the different pieces of equipment at the 
GEOCEIVER STA 31370, a wooden tower was built as shown in 
figure 4.2. This tower was set over the station for the 
period of the test. An NPS equipment vehicle was parked as 
close as possible to the tower. Electrical power was 
obtained from the Monterey Coast Guard facility. 

Prior to the ship test, both GPS antennas were 
placed on the tower and centered over the station mark. For 
the first two nights, 9 and 10 May, equipment problems 
allowed the use of only one system each night and data were 
acquired only to determine the stability of each system. 
During the next night, 11 May, position outputs for the two 
simultaneously operating GPS receivers were recorded. This 
static test gave a determination of any difference between 
the two receivers under almost identical conditions. 

For the ship test, the GPS antenna was centered over 
the mark and was anchored at the top of the tower; the 
geoceiver antenna was centered over the mark and anchored at 
the middle level; and the theodolite on the tripod was set 
over the mark at ground level. This arrangement provided 
colocation of all three systems. It was assumed that no 
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Antenna Tower on Coast Guard Pier 



error was introduced by the presence of the wooden 
structure. The shorebased ARTEMIS antenna unit was offset 
from the theodolite to eliminate interference in the signal 
due to the tower. An offset position was computed relative 
to the station mark and was used in postprocessing of the 
ARTEMIS data. 

Except for the antennas for the GPS, the ARTEMIS, 
and the geoceiver systems, all other equipment was installed 
in the equipment vehicle, where an operator was stationed 
throughout the test to monitor the data quality. 

The GPS Manpack is a manually operated system. As 
mentioned previously, it was decided that manual operation 
of GPS was detrimental to the goal of the test, and that an 
interface that would provide an automated printout of 
positions from the Manpack as well as a printout of GPS time 
would be desirable. Knute Serstis and Gary Whitsell of the 
Office of Testing and Evaluation of Systems (OTES) , NOA A, 
assembled two interfaces that provided data storage on 
cassette tapes as well as the desired printout. Each system 
operated in conjunction with a HP-9825 computer and a 
HP-7845 printer/plotter. Software was provided by OTES, 
NOAA. This arrangement permitted storage of the satellite 
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messages for further processing and provided an automated 
printout of the ARTEMIS data. A GPS Manpack had been sent 
to OTES for testing during the construction of the 
interface. Technical expertise and operation of the 
interface during the test were provided by OTES. 

The data from the GPS receiver were printed at an 
interval of 24 seconds, and consecutive data were 
continually recorded on tape and printed during the test 
periods. In addition, a manual record of the GPS position 
from the CDD was logged as the ' marks' were called via the 
radios. 

The ship positioning systems were installed as shown 
in figures 4.3. and 4.4. An equipment shelter was secured 
to the boat deck of the R/7 ACANIA. A length of six-inch 
diameter polyvinyl chloride (PVC) pipe was attached to the 
aft wall of the building on the centerline of the ship. The 
GPS antenna and the system of lights previously mentioned 
were set on top of the pipe. 

The ARTEMIS mobile antenna unit was set on a 
platform located on the centerline of the ship and at the 
top of the building high enough to not interfere with the 
signal from the ship radar beam. For the period of the 
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R/V ACANIA Equipment Layout 
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Fig. 4.4. Equipment Installation 
on the R/V ACANIA 
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test, the weather conditions allowed ship operations without 
use of the radar, thus there was no possibility of damage to 
the ARTEMIS unit. 

All other equipment — the ARTEMIS receiver unit, the 
GPS receiver, control data unit, power filter, and 
interface, the HP-9825, the HP-7845 plotter/ printer, and 
the real time clock — was placed inside the building. An 
operator monitored the equipment during the ship tests and 
another individual radioed the ’mark', or fix, to the 
theodolite operators on shore. The equipment used at the 
shore location was installed in a vehicle and consisted of 
the same components as on the ship. Operation of equipment 
during the snip tests was monitored by an individual located 
in the vehicle. 

C. LOG OF EVENTS 

The test was divided into two portions — 1) a shorebased 
test which consisted of two nights of stability testing, and 
one night of receiver comparison tests, and 2) the ship test 
consisting of three nights of shipboard operation. The R/V 
ACANIA from the Naval Postgraduate School was the platform 
used for the shipboard tests. 
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1 . Stability Test 



Operations began on 9 May, 1981, GMT. One of the 
GPS systems was placed over the station mark, and continuous 
data was recorded during the satellite availability periods. 
In this manner, system 1 was tested on 9 May, and system 2 
was tested on 10 May. A position was defined as a quality 
observation if four or more satellites were used by the 
receiver to obtain the output. 

2 . Static Test 

Both GPS antennas were centered over GEOCEIVER STA 
31370, and data from simultaneous operation of both 
receivers were recorded at the equipment vehicle. One 
antenna was placed on top of the wooden structure and the 
other immediately below it on the middle level. The 
vertical displacement was assumed to cause negligible error 
in the horizontal positions. Data were acquired on 11 May 
from both receivers with a quality observation defined as a 
position determined from four or more satellites. 

3. Shipboard Test 

For the entire period of shipboard testing, the wind 
speeds were less than 5 knots and the seas were calm with 1- 
2 feet or less of swell. Visibility was excellent throughout 
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the test week. The ship track for each set is shown in 
figure 4.5. A log of subsequent operations follows: 

12 May The GPS Manpack receiver was installed on the R/V 

ACANIA and GPS and theodolite operations were conducted. 
The ship was allowed to drift during sets 1-3 (fix 1-69) 
with 23 fixes per set, in a east-northeasterly direction. 

13 May The ship cruised at 4 knots on track lines from 

south to north for sets 4 and 6, and lines from north to 
south for sets 5 and 7. During fixes 139 through 157 the 
radio at the Coast Guard pier was disabled, and no data were 
acquired by either the theodolite observer or the operator 
in the equipment vehicle. 

14 May The shipboard GPS system was inoperable. The 

trouble was traced to a faulty antenna cable, requiring 
replacement. No data were acquired. 

15 May The ARTEMIS, the GPS and the theodolite network 

were all used for this test period. Data from sets 8 
through 13 were acquired (fixes 165-262). 

16 May The control data unit on the shore based GPS system 

failed before the start of ship operation and could not be 
repaired without factory maintenance. 
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Survey Area and Location of Sets By Day 



A total of 252 fixes were recorded by the theodolite 
network and all GPS positions corresponding to the 
theodolite fixes were manually logged. The ARTEMIS azimuth 
and distance information and consecutive GPS positions on 
the ship and in the vehicle were stored on cassette tapes 
and printed by the HP-7845 printer/plotter. 
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7. DATA PROCESSING METHODS 



Processing of the thesis data was divided into four 
interrelated, but separate, tasks. These tasks were the 
basic geodetic survey, theodolite and ARTEMIS ship position 
determination, differential GPS ship position determination, 
and the statistical analysis of the ship position 
comparisons. Each task will be discussed separately. 

A. BASIC SURVEY COMPUTATIONS 

As described earlier in the Test Procedures, a basic 
geodetic survey of the test area was conducted to place the 
theodolite control stations on the World Geodetic System of 
1972, or WGS-72. To begin the survey, a Magnavox AN/PRR 14 
Geoceiver was placed over the center survey disk, GEOCEIVER 
STA 31370 1981, and the digital Doppler data, from TRANSIT 
satellite 68, were recorded onto paper tape. Thirty-three 
passes, over a period of two and one half weeks, were 
recorded, but analysis of the equipment maintenance 
diagnostics indicated a possible error in the records. 
These data were therefore discarded. 

A second set of forty satellite passes of TRANSIT 
satellite 68 were recorded over a ten day period and were 
determined to be of acceptable quality according to the 
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standards of the DMA Department of Geodesy. This second set 
of data was then submitted to DMA, Department of Geodesy, 
Computation Section, for analysis. Also submitted were 
meteorological data taken during the passes and additional 
information required for documentation of a first order 
TRANSIT position. Satellite 68 Doppler data for the survey 
was then processed in conjunction with the tracked satellite 
position by DMA computer facilities. A final geodetic 
position on the HGS-72 ellipsoid was obtained for GEOCEIVER 
STA 31370. Coordinates of the station are listed in the 
Test Procedures section. The computed height relative to 
the ellipsoidal surface is -33.39 meters. Computation of 
the geographic position utilized a minimum of 38 individual 
TRANSIT satellite passes recorded for GEOCEIVER STA 31370. 
This provided for an accuracy of one part in 10^ (first 
order accuracy standard). [Ref. 38] 

To orient the positioning network, an initial azimuth 
from GEOCEIVER STA 31370 to station MUSSEL 1932 was computed 
from two sets of observations of POLARIS. All observations 
and computations of the star sighting utilized procedures 
outlined in the National Ocean Survey Photogrammetric 
Instructions Number 4, Revision 1, 4/24/73, the Hydrographic 
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Manual, and the Manual of Geodetic Triangulation, SP-247. 
[Ref. 39 , 403 

Having the initial azimuth, distances were measured 
between all three control stations using a Tellurometer MRA 
5 microwave distance measuring unit. Twenty observations of 
distance were made for each of the three lines measured: 
GEOCEIVBR STA 31370 to MUSSEL , GEOCEIVEH STA 31370 to USE 
MON, and MUSSEL to USE MON. The average of the distance for 
each line was then corrected for instrument error and 
meteorological conditions following the instructions 
contained in the Operators Manual and the Tellurometer 
Manual. [Ref. 41 ] Resulting slope distances were then 
reduced to geodetic distances using the respective 
corrections for slope to horizontal distance, reduction of 
elevation to sea level, and the chord to arc distance. 
These corrections were computed using the RGS 1972 
ellipsoidal parameters and the methods outlined in the 
Hydrographic Manual, the Manual of Geodetic Triangulation, 
and Introduction to Surveying. [Ref. 42] The resulting 
geodetic distances meet third- order, class I specifications 
for traverse accuracy, 1:20,000. In addition, angles to the 
other stations were observed at each control station, giving 
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a total of three sets of angles. Each set was comprised of 
four direct and four reverse plate observations with a Wild 
T-2 theodolite. Each pair of direct and reverse 
observations was corrected for initial pointing errors. 
Observations and computations were made using the methods 
detailed in the Hydrographic Manual, the Manual of Geodetic 
Triangulation, and the Manual of Second and Third-Order 
Triangulation and Traverse. [Ref. 43] The resulting 

closure meets the specifications for traverse accuracy of 
three seconds per traverse station or a closure of 6 
seconds. 

With the geodetic position of GEOCEIVER STA 31370 known, 
as well as the azimuth from GEOCEIVER STA 31370 to MUSSEL 
and the distance from GEOCEIVER STA 31370 to MUSSEL, the 
geographic position of MUSSEL was computed on the WGS 1972 
ellipsoid using the direct geodetic computation. This 
computation used the Rainsford's method modified by T. 
Vincenty and including the Helmert's elliptical terms. 
[Ref. 44] 

Next, using the azimuth from MUSSEL to USE MON, the 
geodetic distance, and the results of the direct geodetic 
computations, the geographic position of USE MON was 
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determined. Also, the ARTEMIS geographic position was 
calculated as an offset using an azimuth and distance from 
GEOCSIVER STA 31370. With these computations, the geodetic 
survey was completed. 

Positioning of the ship, R/V ACANIA, required one 
azimuth, or line- of-position, from each of the control 
stations MUSSEL, GEOCEIVER STA 31370, and USE MON. These 
lines of position were used to triangulate the vessel 
position. Angles from each of the three theodolite 
locations were observed as described in the Test Procedures. 
Each of the three angular measurements were then corrected 
for mean pointing errors and the individual theodolite 
collimation errors. Each corrected angle was then adjusted, 
by control station, to obtain an azimuth from south, as 
required to triangulate the vessel position. The three 
resulting azimuths and the corresponding control station 
latitude and longitude were then used in a least squares 
three line of position intersection computation. Results of 
the computation were a vessel position in latitude and 

longitude, error ellipse semi-major radius, error ellipse 
semi-minor radius, azimuth from south of semi-major radius, 
and circular standard error. [Ref. *15 ] Figure 5.1 

illustrates these error ellipse parameters. 
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THE ERROR ELLIPSE 



The mean displacement of a point to be expected in the direction of 
the two coordinate axes, X and Y, due to measuring errors, can be ex- 
pressed by the quantities cr and cr 

* y 

When the north-south axes are rotated, two axes (X and Y) perpendicular 
to one another are obtained, corresponding to the maximum and minimum 
values of cr and cr v . The maximum and minimum values of o’ , and the 
clockwise rotation at^gle (°<) of the major axes of the error ellipse are 
included as part of the output of the Geodetic Package and can be used 
as input if desired. 

The circular standard error (CSE) is also included in the output with the 
error ellipse parameters. The value of CES is the radius of a circle 
whose area is equal to the area of the error ellipse. That is: 



CSE 




tf" 

min 




Fig. 5.1. 
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For this study, 262 vessel positions, or fires, were 
observed and recorded. On examination during processing, 51 
fixes were found to contain errors and data omissions which 
rendered the vessel positioning computation invalid or 
indeterminate. The remaining 211 fixes comprised the data 
set which was then broken down into smaller sets for 
analysis and comparison. These sets were numbered 1 through 
13 and categorized according to vessel track. All sets were 
then processed using several different methods to determine 
the most accurate approach. 

3. THEODOLITE AND ARTEMIS SHIP POSITION DETERMINATION 

Before processing the GPS Manpack data, the ARTEMIS data 
were processed to ensure the vessel positions, as obtained 
by theodolite, were correlated with the correct automated 
and hand recorded GPS data. ARTEMIS data consisted of an 
angle, MUSSEL-ARTEMIS OFFSET- vessel, and distance from the 

'•v 

ARTEMIS OFFSET to the vessel. The angular measurement was 
first converted to an azimuth from south using the azimuth 
of ARTEMIS OFFSET to MUSSEL as the correction and adding the 
measured angle. Vessel position based on ARTEMIS was then 
computed using the direct geodetic program (also used for 
the basic survey) . Comparisons of the resulting geographic 
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positions and the theodolite derived geographic positions, 
were made for corresponding fix numbers using the inverse 
geodetic computation on the WGS-72 ellipsoid. The inverse 
computation utilizes two geographic positions as inputs for 
calculating the forward and back azimuths and the distance 
between the positions. [Ref. 46] For comparison, if the 
two positions differed by more than 20 meters, the ARTEMIS 
positions on either side of the suspect ARTEMIS fix were 
then used to compute the inverse computation to the 
theodolite ship position. In seven fixes, the fix position 
values were changed in the May 15 GPS data to give more 
accurate results. 

C. DIFFERENTIAL GPS SHIP POSITION DETERMINATION 

With the check and adjustment of the GPS data and the 
theodolite determined vessel position complete, the first 
GPS data processing was the comparison of the uncorrected, 
non-differential ship GPS Manpack receiver data to the 
theodolite computed ship position. Using the inverse 
geodetic computation, the forward azimuth and the distance 
between the theodolite computed ship position and the 
ncn-dif ferential GPS ship position were completed. These 
results were then used in the analysis of the differential 
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processing to examine the stability of the shipboard Manpack 
receiver compared to the GPS Manpack receiver located over 
GE0CEI7ER STA 31370. 

There were two methods of analysis of the differential 
GPS. The first of these was the determination of a static 
correction to be applied to all the GPS data. This static 
correction in the differential GPS computations was to 
account for a suspected bias between the two similar, but 
not identical, Manpack receivers. A mean static correction 
was computed for the most stable reception period during the 
static tests on 11 May. With the two receivers colocated 
over the GE0CEI7ER STA 31370, a mean correction of -0.954 
seconds of latitude and -0.307 seconds of longitude was 
found between GPS Manpack System 1 and System 2 during times 
07:16:00 GMT and 07:25:12 GMT. This static correction was 
then added to the latitude and longitude of GPS Manpack 
System 1 (located over GE0CEI7 ER STA 31370) data. 

The second method of analysis or processing in the 
differential mode used the latitude and longitude obtained 
directly from the GPS Manpack receivers. la either 
processing method, the computation of the differential 
corrector to the shipboard GPS Manpack receiver data was the 
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same. 



An explanation of the computations and processing of 



the differential correction for both methods follows. 

Computation of the differential corrector for each GPS 
fix position was accomplished by using the geographic 
position obtained with the stationary Hanpack receiver, 
located over GEOCEIVER STA 31370, and the geographic 
position of GEOCEIVER STA 31370. These two positions were 
then used in the inverse geodetic computation to obtain a 
differential correction vector (forward azimuth, and 

distance) from the stationary GPS receiver to GEOCEIVER STA 
31370. In figure 5.2, the differential correction vector is 
labelled "A 1 ’. Osing the direct geodetic computation, the 
differential correction vector was applied to the 

corresponding shipboard GPS Hanpack receiver position, and 
the differential GPS position was obtained. These resulting 
positions were computed on a one-to-one basis using the fix 

/N 

numbers as identifiers. Vector ’'B”, the differential 

corrector vector which is identical to the vector "A" in 
figure 5.2, is applied to the shipboard GPS Manpack receiver 
and results in the differential GPS ship position. 

Comparison of the differential GPS ship position with 
the theodolite-computed ship position was accomplished using 
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the inverse geodetic computation for both methods of 
differential GPS processing. The resulting forward azimuth 
and distances for both processing methods were used in the 
statistical analysis. 

D. STATISTICAL DATA ANALYSIS 

Statistical analysis was done on data separated into 
sets according to the method of data processing. The static 
data and stability data were graphed as functions of 
latitude and longitude versus time. User range error, 
furnished by the Space Division, GPS/OL-AQ, Vandenburg, AFB, 
was plotted as a function of time for each day of the test 
(see Results) . Other graphed relationships are the 
gecgraphic positions of the stationary GPS Manpack receiver 
data versus the GEOCEIVER STA 31370 and the theodolite 
semi- major radius as a function of fix number (see Results) . 

Other statistical data analyses were represented by 
histograms. Plotted on these histograms are the relative 
frequency on the left margin, variable bar increments scaled 
on the bottom, number of variable observations in each bar 
increment on the top, and the total sample size on the top 
margin. Within the histogram the frequency bars are denoted 
by asterisks, the mean by a capital M , the quartiles by 
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periods and the empirical density function by the capital F 
curve. The empirical density function was computed using: 



F(Z) 



1 N (X f Z) 

£ JJ 1— 

N (B (N) ) i= 1 B (N) 



[ Bef . 47] where: F(z) = the empirical density function, z 

= the frequency of observation, x = individual observation, 
N = the total number of observations, H = weight function, B 
= bandwidth, 

W = 1 - z 

B (N) = range of all observations / <7 
The statistical values listed below the histogram are 
defined in Table V. Computed histograms show the distances 
between: 1) the theodolite-computed ship position and the 

differential GPS ship position, both with and without the 
static correction, by set; 2) the theodolite-computed ship 
position and the differential GPS ship position, without 
static correction, for five satellite observation periods on 
each of the three ship test days; 3) the theodolite-computed 
ship position and the non-differential GPS ship position for 
the five satellite-observation periods on each of the three 
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TABLE V. STATISTICAL DEFINITIONS 



The definitions of the basic statistics computed by the program are 
listed below:. 

MEAN: Average of the sample 

MEDIAN: Mid-Value of the sample, if there are an odd number of 

sample points, or the average of the two middle values 
for an even number of points 

TRIMEAN: 0.25 * (C^ + 2Qg + Qg), where are the quartiles 

MIDMEAN: The average of all sample values lying between the upper 

and lower quartiles 

MIDRANGE: Average of the maximum and minimum 

GEOMETRIC MEAN: GM = n *Jx 1 -x„*. . . x 

1 n2 n 

HARMONIC MEAN: HM = n / ^ 1 /x i 

VARIANCE: Unbiased estimators for variance and standard deviation 

are used - th e square of the standard deviation 

STANDARD DEVIATION: S = \.$? x? - n x 2 - 

V - 1 

COEFFICIENT OF VARIATION: Standard deviation / | Mean | when the 

mean is less than 10. ** -30, the coeffic- 
ient of variation is set to zero. 

MEAN DEVIATION: The average of the sum of the absolute differences 

between the sample values and the median 

RANGE: Maximum - minimum 

MIDSPREAD: - Q-j , also called the interquartile distance 

M3: Third central moment. Unbiased estimator is used 

M4: Fourth central moment. Unbiased estimator is used 

COEFFICIENT OF SKEWNESS: M3 / (Standard Deviation) ** 3 

COEFFICIENT OF KURTOSIS: (M4 / (Standard Deviation) ** 4) - 3 

BETA1 : Biased estimate of third central moment; can be used in 

testing for normality 

BETA2: Biased estimate of fourth central moment 

MAXIMUM: Largest sample value 

MINIMUM: Smallest sample value 

QUANTILES: The alpha-quantile, x(alpha), is the solution to the 

equation: 

Probability (x .le. x(alpha) ) = alpha 
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ship test days; and 4) the GEOCEIVER STA 31370 position and 
the uncorrected stationary GPS position for the three days 
of ship tests. In addition, the azimuths between the 
theodolite-computed static corrections, for each subset, 
were graphed as a series of histograms. To illustrate the 
difference between the GEOCEIVER STA 31370 position and the 
stationary GPS position, the azimuths between the two 
positions are graphed on histograms for each of the three 
test days. 
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VI. ERROR ANALYSIS 



The error sources contributed by each positioning system 
are described. Estimates of the effect of each on the final 
results are determined where possible. 

A. ARTEHIS POSITIONING SYSTEM 

The sources of error in the ARTEMIS ship position can be 
divided into two categories — those sources of error inherent 
in all microwave positioning systems and those unique to the 
ARTEMIS system. All are investigated and the contribution 
of each to the total error is determined. 

1 . Micro wave System Error s 

Sources of system error include range holes, 
inaccuracies in calibration, non-ideal geometry of the lines 
of position describing the fix, and interference with 
reflections of the signal from nearby objects. 

Range holes are a cause of unreliability in 
microwave positioning systems. The problem exists when the 
sea surface is calm enough to reflect part of the 
transmitted signal. The reflected signal will arrive at the 
receiving antenna at a later time than the signal which 
travels directly, causing a phase difference between the two 
signals. A zone or region of destructive interference of 
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the directly transmitted signal and the reflected signal 
will occur within which the vessel will receive an 
attenuated signal. When the phase difference of the two 
paths is exactly 180 °, no signal will be received provided 
the reflected signal is as strong as the direct-path signal. 
Prediction of the range holes, therefore, is useful in 
survey planning so as to avoid this problem. 

The range hole is a function of shore unit and ship 
unit antenna heights, the wavelength of the system, and the 
distance from the shorebased unit. If the operating area is 
known, antenna heights can be determined and equipment can 
be installed so that no range holes exist in the survey 
area. If fixed antenna heights are known, the region of 
uncertainty can be determined from the equation: [Sef. 48] 

The following equation is used for this computation: 

D n = 2 (hi ) (h2)/nX, 

where X is the wavelength (constant for ARTEHIS at .032 m) , 
hi and h2 are antenna heights, and D n is the distance from 
the shore based unit to the nth range hole. For this study, 
the antenna height of the fixed unit above mean sea level at 
the GEOCEIVER STA 31370 was 3.8 meters. The mobile antenna 
was installed on top of the equipment shelter on the boat 
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deck of the R/V ACANIA at a height of 7.9 meters above the 
the waterline. Two range holes within the survey area were 
determined to be located at ranges of 1876 meters for n = 1 
and 938 meters at n = 2 (fig. 6.1). These two ranges were 
computed for the case of mean lower low tide. The change in 
width of the range hole due to a tidal fluctuation of 1.5 
meters during the test is shown in figure 6.1. (During the 
night of 15 May the sea surface was sufficiently disturbed 
by a two-foot swell and a surface wind of 2-4 knots such 
that surface reflection was minimal and no loss of signal 
was observed at the computed ranges) . 

Calibration of the ARTEMIS system was accomplished 
in two stages. As part of the study, the ARTEMIS was 
obtained for a one week test of the system from 13 to 17 
April. During this time it was calibrated for azimuth and 
range over the baseline between the GEOCEIVER STA 31370 and 
MOSSEL. The azimuth readings required no adjustment but the 
range was reduced by two meters. The equipment was returned 
to Marina v Corporation to fulfill prior commitments and then 
was returned to Monterey for the ship test on 15 May. 
Before the evening operation, the ARTEMIS was used to 
measure the same baseline as previously. A plus two meter 
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Fig. 6.1. ARTEMIS Range Holes 



bias was observed at this time, but no immediate 
compensation was made. This error was assumed constant 
throughout the test and was included in the calculation of 
the ship-determined ARTEMIS position during postprocessing. 

One of the advantages of the ARTEMIS system, unlike 
other microwave positioning systems, is that the geometry is 
always ideal for the lines of position that determine the 
fix. The range-azimuth technique of determining the vessel 
position defines an uncertainty ellipse having axes that are 
always mutually perpendicular. This allows prediction of 
errors as a function of range only, independent of the 
angle. 

The manufacturer's specified error in range 
measurement of the ARTEMIS is + 1.5 meters for all line of 
sight ranges. One axis of the uncertainty ellipse is 
therefore constant. The other axis, that describes the 
error introduced by the system azimuth uncertainty, is 60 cm 
per 1000 m range. (This is determined from an angular error 
of 0.033 ° or 2' of arc) . Both system error values represent 
an accuracy of twice the standard deviation value, or a 95^ 
probability of the true value occurring within the margin 
allowed by these parameters. 
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The positional uncertainty can be expressed in terms 
of the radius of a circle of equivalent probability within 
which the true fix is expected to occur. The method of 
computation of this radius is described in Appendix A. The 
technique allows the conversion of uncertainty ellipses to 
radii of circles of equivalent probability. Values from 
this computation for the radii of circles of 90S probability 
at different ranges used in the survey area for 15 May are 
displayed in figure 6.2 along with the circles of 90S 
probability, the magnitude of the radii at the ranges used 
for the study do not exceed 5.4 meters. One goal of this 
study was to measure the difference in positional accuracy 
between ARTEMIS and GPS. A maximum error of 1.0 meters, an 
order of magnitude smaller that the tested system accuracy, 
is desirable. At the present time, however, there are few 
systems that can provide the order of accuracy system for 
the position of a moving vessel, so the 5.4 meter error 
observed at 4000 meters had to suffice. 

Error caused by reflection of the signal from nearby 
objects is difficult to measure. To minimize this problem, 
the offset position for the ARTEMIS unit was planned so as 
to be clear of possible interference from the wooden 
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structure and the ejuipment vehicle. No signal interference 
was observed despite the ship's motion. This source of data 
unreliability was therefore not considered to be 
significant. 

2. ARTEMIS System Errors 

Errors unique to the ARTEMIS system include a 
pointing error source and the lack of resolution of the 
azimuth measurement in the dynamic mode of operation. 



The operation of the ARTEMIS system requires a 



manual alignment of 


the 


fixed 


antenna 


SO 


that it 


is 


perpendicular 


to a known 


object. 


This 


is 


necessary 


to 


establish an 


original 


reference 


direction 


upon which 


all 



angular measurements are based. It is possible that a 

pointing error could be produced by the operator if the 
crosshairs of the telescope used to align the antenna are 
not exactly sighted. This would introduce a systematic 
error into all subsequent angle measurements. To minimize 
this source of error, the alignment of the antenna was 
checked before and after the period of operation. The 
original azimuth displayed on the antenna unit display 
agreed exactly with the final azimuth value. Any 
contribution to ship position error was considered minimal. 



no 



The ARTEMIS was operated using the dynamic mode 
option (see Equipment Descript ion--ARTEMIS) . The resolution 
of the system in this mode is 1 meter for displayed range 
and 0.1° for azimuth values. This resolution limit 

i O 

introduces an uncertainty of _0.5 meters in range and 0.05 
in angle measurements. The computed uncertainty ellipse as 
expressed by radii of circles of equivalent probability, is 
shown in figure 6.2. The system error and all other errors 
mentioned above are assumed to occur within the scope of 
this resolution uncertainty. 

Another source of system error has been proposed in 
a study conducted in June, 1980, by NOAA. [Ref. 49] The 
test compared the ARTEMIS system with two ranges measured by 
two Del Norte Trisponders (a short range microwave 
range/range positioning system) and two lines of direction 
measured by two shore-based theodolites. The results 
described an underway bias of -0.05° of unknown origin in 
the angle measurements. This discrepancy has not been 
resolved. The preliminary test conducted in Monterey in 
April, 1981 compared the ARTEMIS system with three lines of 
direction determined by three theodolites. The purpose of 
the test was to resolve the underway bias problem. Results 
are enclosed in a separate report. [Ref. 50] 
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B. CONTROL NETWORK 



The entire survey was oriented relative to the TRANSIT 
position on the Coast Guard Pier (GEOCEIVER STA 31370). 
Therefore, errors inherent in the position of this station 
mark affected all results. Assuming that the geographic 
coordinates published by DMA/HTC represented the best 
estimate available of the GEOCEIVER STA 31370, the remaining 
significant sources of error affecting all results were 
survey errors. The control network for the study was 
determined from the traverse and it was assumed that all 
positions used in the study were oriented relative to each 
other on the WGS-72 ellipsoid. A change in the ellipsoidal 
surface within the area encompassed by the traverse was 
assumed to be insignificant. 

Survey errors stem from the methods used for obtaining 
the initial azimuth for orientation of the survey, the type 
of traverse chosen, and the geometry of the traversed 
triangle. 

To obtain a starting azimuth for the traverse, two sets 
of observations on Polaris were recorded at the GEOCEIVER 
STA 31370. A striding level attachment to the theodolite 
was not used. The instrument level bubble was observed 
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repeatedly during the observations to ensure constant 
instrument alignment with the vertical. Failure to keep the 
instrument level could have introduced a + 30” uncertainty 

into the starting azimuth of the survey. [Hef. 51] To 
consider the significance of such an error, all stations 
were recomputed using a deviation of + 15” from the computed 
azimuth. The resulting shift in the positions of MUSSEL and 
USE MON are listed: 

USE MON 0.133 m distance at 230° 36 ' 47 . 16” azimuth 

MUSSEL 0.132 m distance at 230° 33' 4.08” azimuth 

The change is considered insignificant relative to the other 
errors described above. 

The survey was a closed loop traverse using SE0CEIVER 
STA 31370 as the origin and closing position. The weakness 
in this method involves the measured distances. A method of 
checking the observed angles in the triangle exists via the 
angular closure if a significant deviation in one of the 
angles can be detected. The distance measurement, however, 
can produce an undetected bias in the survey if the distance 
measuring instrument contains a systematic error. This 
error would shift all the surveyed positions. Correct 
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measurement of distance by the MRA V Teliurometer , 
therefore, was important. The equipment was calibrated over 
a known baseline and the errors were within the centimeter 
range. This error was therefore negligible. 

The geometry of the triangle was not optimum. Ideally, 
the triangle would have been equilateral to minimize the 
source of observing error. The choice of traversed points 
was based on the location of preexisting station marks and 
the desired geometry of the theodolite fixes to the ship. 
The closures computed for the traverse met third order 
specifications for accuracy, and the survey was not adjusted 
to remove the closure error. [Ref. 52, 53J 

The survey errors were considered to be within the 
uncertainty cf the TRANSIT station location and were 
therefore not considered significant for the study. 

C. TBEODOLITE NETWORK 

The errors generated by the theodolite network are 
divided into equipment errors, operator related errors and 
system errors, of which the latter is the most significant. 

The equipment used for the theodolite network were three 
1-second theodolites which have a measurement resolution of 
10.5". All theodolites were in good condition. 
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Instrumental errors of a theodolite are listed in Table 71. 



[Bef. 54] 

The major source of equipment error was introduced by 
not changing face, the process of observing angles a second 
time with the telescope reversed. Due to the nature of the 
test coupled with 1-minute fix intervals, it was not 
possible to observe all the fixes in both the direct and 
reversed modes. As an example of the magnitude of this 
error, the horizontal collimation errors were determined for 
each theodolite prior to the study. The results are listed 
below: 

Theodolite Collimation 5rror 

# 14405 - 1.0“ 

# 14482 + 2.7” 

# 14452 - 0.7” 

All the measured angles were adjusted for this error as 
described in the processing section. The largest influence 
on the remaining instrumental errors was the non-vert icality 
of the main axis. If the non-verticality were large (for 
example, one division off-level in the lower spirit level 
bubble), the angle measured could be in error by 7”. 
[Bef. 55] To minimize this error, each theodolite was 
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TABLE VI. Theodolite Instrument Errors 



ERROR 

Eccentricity, misplaced indices 

Graduation of circles 
Micrometer run and vernier error 

Non-vertical i ty of main axis 

Bad setting of the upper spirit 
1 evel 

Vert i cal ity of cross wires 

Horizontal collimation 

Vertical collimation and index 
setting error 



ELIMINATED BY: 



Reading both sides of circle 
and changing face 

Changing zeros 

Selection of suitable zero 
setting 

Error is not cancelled by any 
observing procedure 

Correct setting or applica- 
tion of correction 

Changing face and using same 
part of wire 

Changing face 

Changing face 
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relevelled if necessary between sets of observations. 



The 



magnitude of this and the collimation errors were considered 
to be insignificant relative to the system errors discussed 
later in this section. 

The operator-related errors are a function of the 
individual performance of the observers. The experience and 
personal comfort of the observer, for example, are two major 
factors which can influence the observed angles. An 
observer's performance and resultant impact on the quality 
of the results depend on the judgment and reaction time of 
the observer and are influenced by variables such as fatigue 
and age. The subjective nature of the problem makes an 
estimate of the error difficult, therefore, for this study 
the error was considered random. 

The major source of error in the theodolite positioning 
system dealt with the timing of the fix. A delay existed 
between the time the 'mark' was called over the radio and 
the actual time the angle was observed. The timing problem 
was most critical when the ship was close to one of the 
three observer locations because the rate of change of the 
bearing of the ship was greatest at this time. Hence a 
small delay generated a large amount of error in the angle 
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recorded. To indicate the effect this error source may have 
had on the position of the ship, two fixes were chosen 
arbitrarily — one at close range to each of MUSSEL and 
GEOCEIVER STA 31370. A one minute value was added to the 
angle which was most affected by the timing error and the 
fixes were recomputed. The same process was repeated for 
the three fixes with the one minute value subtracted from 
the critical angle in each fix. The results follow: 

Number Original Fix * 1 Minute - 1 Minute 

210 36°37* 20. 6 22 82" 36° 37* 20.6 146 8" 36 ° 37 » 20 . 63081 " 

121°53« 41. 97184" 121° 53* 4 1 .96933” 1 2 1 ° 53 • 4 1 . 97426" 
259 36°36»45. 88048" 36° 36« 4 5. 87574" 36 °36 * 45. 89082" 

121°52* 54.05062" 1 2 1 ° 52* 54 . 046 35 " 121°52» 54. 05459" 
No significant error was generated from this source. 

Another source of system error was associated with the 
geometry of the fix. Ideally, the three theodolite angles 
would describe two equilateral triangles to the ship; 
however, within the operating area, the geometry varied. As 
an indication of the size of this error, representative 
fixes were chosen and the circular map accuracy standards 
(CMAS) of each of the combinations was determined. Figure 
6.3 illustrates the approximate error caused by variations 
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Fig. 6.3. Representative Theodolite Intersections for Error Ellipse Computations 



in the geometry. Nithin the area outlined in the figure the 
geometry creates a CMAS of not more than 10.4 meters. This 
is therefore the source of the largest probable error in the 
theodolite method of positioning. Nevertheless, this method 
in combination with the ARTEMIS system provided the best 
estimate of the ship location. 

D. TRANSIT POSITIONING SYSTEM 

The method used to determine the TRANSIT position 
(GEOCEIYER STA 31370) is the most accurate satellite 
positioning technique presently available. Error sources 
are reduced both during data collection and during 
postprocessing of the recorded values. 

DMA places requirements on the collection of data that 
will minimize error. For this study, only satellite passes 
whose maximum altitude ranged between 10° and 30 0 were 
recorded. Meteorological data (temperature, pressure, and 
humidity) were collected for each pass, and a determination 
of the receiver clock drift was computed after each pass. 
The amount of clock drift as a function of time gave a 
real-time indication of receiver performance. Equipment 
checks were conducted routinely, and the results were 
included with the data package submitted to DMA. 
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An error budget of the TRANSIT system for point 
positioning is listed in Table VII. [Ref. 56] The precise 
ephemeris method of postprocessing accounts for all of these 
errors: 1) for each pass the tropospheric correction is 
computed from the meteorological data; 2) the ionospheric 
correction is determined by observing the shift in frequency 
between the two transmitted frequencies; 3) the receiver 
clock error is known; 4) the satellite orbits are observed 
by the tracking stations so that precise satellite positions 
are known; 5) errors due to uncertainties in the 
geopotential model and surface drag forces are accounted for 
during computation; and 6) the receiver altitude relative to 
the ellipsoidal surface is computed approximately from the 
given MSL elevation. 

The accuracy of the computed TRANSIT position is stated 
as 1.5 meters per axis within a 90S confidence interval. 
This is assumed to provide the best estimate possible of the 
G EOCEIVER STA 31370 relative to the HGS-72 coordinate 
system. 

E. GPS ERROR SOURCES 

Global Positioning System error sources are separated 
into two main categories, the space vehicle system errors 
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TABLE VI I. TRANSIT System Errors 



SOURCE ERROR (m) 



Uncorrected propagation effects 
(ionospheric and tropospheric effects) 1-5 

Instrumentation and measurement noise 
(local and satellite oscillator phase 
jitter, navigator's clock error) 3-6 

Uncertainties in the geopotential model 
used in generating the orbit 10-20 

Incorrectly modelled surface forces 
(drag and radiation pressure acting on 
the satellites during the extrapolation 
interval) 10-25 

Uncertainties in navigator's altitude 
(results in a bias in longitude) 10 

Ephemeris rounding error 5 
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and the GPS user receiver errors. 



There are seven possible 



sources of error in the space vehicle system--errors 
associated with satellite equipment, satellite clock, 
ephemeris, atmospheric delays, signal multipath effect, 
satellite update, and satellite propulsion firing. Table 
III lists these error sources and the approximate user 
equivalent range error for each. [Bef. 57] 

The satellite equipment error, also termed the group 
delay, is caused by differences and uncertainties in 
processing and transmission of the message through the 
individual satellite's circuitry. Each satellite has the 
group delay calibrated prior to launch and this known group 
delay is included in the satellite time corrector within the 
navigation message. However, unknown delays in signal 
processing and circuitry not accounted for by the corrector 
contributes an estimated one meter to the error budget. 

Satellite clock error is caused by variations in each 
satellite clock. Satellite clock time may depart by up to 
976 microseconds from the correct GPS time. Corrections for 
the predicted clock variations are included in the 
navigation message. Residual clock variations contribute to 
the error budget and are included in ephemeris errors since 
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the clock errors are small and have the same effect as 



ephemeris errors. 

Ephemeris errors are the result of satellite clock 
variations, drift by the monitor station clocks, monitor 
station signal processing delays, lack of precise monitor 
station geodetic positions, solar wind and pressure 
parameter variations, and earth gravity model error. These 
sources all contribute to the ephemeris error budget. The 
ephemeris error budget is monitored over long time periods 
and translated by the Master Control Station into an 
User-Equivalent Range Error (ORE). When the ORE erceeds the 
prescribed value of 4 meters the satellite ephemeris data is 
uploaded to reduce the ORE. [Ref. 58] Most of the ORE is 
included into the user clock corrector obtained in the GPS 
position computation. However, not all of the ephemeris 
error is removed during the position computation, and the 
residual error has been estimated to be approximately 1.5 
meters. See Table III. 

Atmospheric delays in signal propagation are the result 
of two effects, ionospheric delay and tropospheric delay. 
Of these, the ionospheric delay, which depends on the 
frequency of the transmitted signal, is greater. With the 
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two frequencies, Li and L2, transmitted by each satellite to 
the user receiver, the ionospheric propagation delay can be 
modelled. Using this model, the majority of the propagation 
delay is accounted for, but a residual delay results in an 
error. The other part of the atmospheric delay error 
depends on the tropospheric propagation delay. Propagation 
delay in the troposphere is a function of signal path length 
through the layer, the receiver air temperature, and air 
pressure and humidity. Since the tropospheric delay is not 
a function of frequency and is highly localized, the 
tropospheric delay corrector has not been incorporated into 
the atmospheric correction computation and is therefore an 
additional source of error. These two error sources combine 
to yield an atmospheric delay error estimated to have a 
magnitude of between 2.4 to 5.2 meters. 

Ionospheric scintillation of the GPS signal, in addition 
to the ionospheric propagation delay, cause errors in 
positioning. Scintillation has the effect of random 
interference with the propagated signal as the signal is 
transmitted through the upper ionosphere. [Ref. 59] 
Interference caused by scintillation combines with the GPS 
user receiver processing noise to cause a random ranging 
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interruption of 



error and, if the interference is severe, 
the pseudo-random number binary code causes loss of signal. 
Errors and signal losses caused by scintillation cannot be 
predicted or estimated. 

Signal multipath errors result from the reception of a 
satellite* s signal transmitted via different propagation 
paths. Received signals thus may have been distorted or 
reduced in strength to a sufficient degree so as to affect 
the range determination from the satellite. Error due to 
multipathing is highly localized and cannot be modelled. 
Estimates of the error magnitude are between 1.2 and 2.7 
meters. 

The last two space vehicle system error sources are 
attributable to the Master Control Station (MCS) satellite 
accuracy and maintenance functions. Usually the MCS uploads 
each satellite at a given time each day and the ephemeris 
data is changed (Table VIII). Should the user not be aware of 
occurrence of upload and if the receiver therefore does not 
obtain the new ephemeris data, ranging to the satellite in 
question will contain a possibly very large error. In 
addition, if the MCS, through the Remote Stations, 
determines that the four skewed reaction wheels cannot 
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TABLE VIII. EPHEMERIS UPDATE TIMES 
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May 15 0515 0445 0440 0545 0545 

0645 0715 

0815 



maintain satellite antenna alignment with earth. 



the 



hydrazine propulsion system is activated and the satellite 
is rotated. Once the propulsion system has used, the 
satellite orbit cannot stabilize immediately, and the 
ephemeris data will be inaccurate. A period of up to 
several hours may be required for orbit stabilization and 
for the new ephemeris data to be computed and uploaded. 
During that period, the ranging computations to the 
satellite will have large errors. The magnitude of these 
two error sources cannot be estimated but is assumed to be 
several orders of magnitude larger than the previously 
stated errors. 

The second error category is that related to the GPS 
user receivers. These sources of error associated with the 
user GPS receiver, the Manpack, are receiver signal noise, 
receiver vessel velocity, signal resolution, method of 
signal processing, ephemeris update capability, and, for the 
differential GPS computations, system accuracy comparison. 

Received signal noise errors result from interference 
introduced into the receiver GPS signals through processing 
and internal receiver circuitry. Additionally, quality of 
signal resolution is dependent upon receiver type and can be 
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a source of error. 



The combined errors from these two 



sources are estimated to be 1.5 meters. [ Ref . 60 3 

Estimates of these error sources for the HVUE have not been 
made but are estimated to be of the same magnitude. 

Velocity of the receiver vessel results in GPS errors of 
varying magnitude. The higher the velocity, the larger the 
error. This error is due to the time lag required to obtain 
and process each NAVSTAR satellite range with the MVOE 
single channel receiver, to compute the average vessel 
velocity and to use the average velocity for positional 
computations. The Manpack requires a six second time 
interval to obtain and process the pseudo-range from each 
satellite used. Thus the fix position is an extrapolation 
of the four satellite ranges to obtain a convergent fix 
position. At one knot, the vessel travels 0.51 meters in 
one second. During the full twenty-four seconds required to 
obtain four satellite ranges, the vessel would have 
travelled 12.24 meters. With the vessel travelling at four 
knots, the distance travelled during the total time interval 
required for a GPS fix position would be 48.96 meters. 
Processing the pseudo-ranges from four satellites requires 
the solution including the average vessel velocity during 
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the 24-second ranging period. Ose of this average neglects 
the relative velocity of the receiver antenna which is 
mounted on the oast. During this test, the GPS antenna was 
mounted 9.06 meters above the waterline of the R/V ACANIA . 
Assuming five degrees of pitch and roll and a vessel roll 
period of 15 seconds, relative average antenna velocity was 
,0.16 meters per second. The antenna velocity error 

translates into a random error of 0.16 meters. The other 
velocity errors have not been estimated but are assumed to 
be less than one meter for the low vessel velocities used in 
the test. 

One error source mentioned previously requires further 
clarification: error due to the relative position at the 

time of the fix. The Manpack requires up to twenty-four 
seconds to obtain and process the satellite ranges. The 
position of the fix and the associated time of fix are then 
recorded. However, the time of the fix may not have been 
the same time the computational results were printed. Thus 
the output position may not have been correlated with the 
correct theodolite computed ship position since the 
theodolite lines of position were observed when the GPS 
position was printed. Errors from this bias source may 
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range from 12.24 meters at one knot of vessel velocity to 
48.96 meters at four knots. This error varied with the 
processing time required to obtain a convergent position. 
An estimate of this error magnitude was not made because of 
its wide variability. 

The type of receiver and the method of signal processing 
also contributed to the error. Hultichannel receivers 
designed for Phase I testing have a lower magnitude error 
for satellite signal reception and processing than do the 
single channel receivers, such as the Hanpack. In the 
multichannel receivers, four or more satellite pseudo-range 
data are taken simultaneously. The single channel receiver 
sequentially receives each satellite and computes the 
pseudo-range data. Simultaneous data reception reduces the 
user velocity error and therefore has less position aliasing 
than that found in the single channel receivers. 

The method of signal processing may induce positional 
errors through the use of non-optimal weighting functions in 
the position computation. Should the ephemeris or ranging 
data change significantly, as happens during uploading, 
signal processing using Kalman filters as used in the 
aanpack do not change the computed position rapidly enough 
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to reflect the large changes in ranges and result in a slow 
change drift of the GPS position. Estimated time 
requirements for user receivers to stabilize after an upload 
are up to two hours. [Ref. 61] 

The capability of the user GPS receiver to receive 
ephemeris data immediately after upload reduces the 
ephemeris errors caused by using an outdated ephemeris. The 
ephemeris data is repeated every thirty seconds and its 
length requires sections of two message subframes. 
[Ref. 62, 63H. The use of current ephemeris data in 
the computations reduces the error for the ephemeris to the 
value shown in table III. The Man pack has the capability to 
accept new ephemeris data immediately, but during the test 
large errors were introduced into the computation on 
satellite ephemeris upload because this selected optional 
function either did not operate properly or was neglected by 
the operators. 

For the differential GPS user, an additional error 
source is introduced into the error budget: the assumption 
that both GPS receivers will observe, obtain the 
pseudo- ranges , and compute the position solution in exactly 
the same manner. If this were true then both receivers 
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should be able to compute the identical position if placed 
at the same position at the same time. With the Manpacks, 
the assumption that the units were identical and the 
position solutions were identical was not made because prior 
testing established that the receivers differed in position 
computations by as much as 3.08 meters. [Ref. 64] 

For the static testing on 11 May, one error source not 
previously discussed was the vertical separation of the two 
Manpack antennas, colocated over the GEOCEIVER STA 31370 
control station, and the possible interference of the 
uppermost antenna. System 1, with the GPS signals received 
by the lower System 2 antenna. Separation of the two 
antennas was 2.1 meters. Calculations using the speed of 
propagation in a vacuum and the two GPS frequencies produce 
wave lengths of 0.19 meters and 0.24 meters. Since 
separation of the antennas was much greater than either 
wavelength theoretically there should have been no 
positional errors due to antenna interaction or 
interference. Differing elevations of the antennas above 
the ellipsoidal surface also may have induced errors in the 
position computation. Estimated errors for the elevation 
differences are less than 0.1 meters. Total errors for the 
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antenna vertical elevation separation and antenna 
interference or interaction have not been estimated. 

Thus, the error budget for the Global Positioning System 
suggests a combined error of up to 57.35 meters. (Jsing the 
differential technique, the bias errors can be reduced in 
magnitude or eliminated. These bias errors include 
ionospheric propagation delays, tropospheric propagation 
delays, variations from predicted ephemerides, and satellite 
clock perturbations. Teasley, Hoover, and Johnson (1980), 
evaluated the error budget for a more advanced GPS receiver 
and found the predicted filtered solution error was reduced 
from approximately 11 meters to approximately 3 meters with 
differential GPS processing. [Ref. 65] Error budgets for 
the tianpack and the differential correction vector 
processing method used in this test did not attain the four 
meter accuracy but were not expected to because the 
receivers used were less accurate, and of different types. 
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VII. RESULTS 



Results of the GPS tests have been separated into four 
sections: stationary GPS evaluation, evaluation of non- 
differential GPS, differential GPS evaluation with static 
correction, and evaluation of differential GPS without 
static correction. Tables IV and VIII list the optimum 
satellite observation periods and the NAVSTAR upload times 
for each day of the testing period. 

A. STATIONARY GPS EVALUATION 

Stationary GPS receiver evaluation was separated into 
the stability tests, the static tests and the tests of the 
stationary GPS receiver used during differential testing. 
Figures 7.1, 7.2, and 7.3 illustrate the results for the 
stability test conducted on 9 and 10 May. The receiver was 
less stable than expected and depended to a great extent on 
the NAVSTAR uploading times and the URE of NAVSTAR 1 (fig. 
7.4 and 7.5). Correlation between times for the NAVSTAR 
upload and changes in the stability can be observed at 0620 
on 9 Say (fig. 7.1), 0655 on 9 Say (fig. 7.2), 0815 on 10 
Say (fig. 7.3), and 0845 on 10 May (fig. 7.3). The 
relationship between stability and URE is illustrated by 
comparison of the URE for NAVSTAR 1 at 0640 on 9 May (fig. 
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GREENWICH MEAN TIME (HOURS) 



7.4) with the stability plot on figure 7.1 for the same time 
period. On figure 7.2 the stability plot variation at 0745 
on 10 May remains unexplained but may have been due to a 
receiver malfunction. 

Results for tests on 11 May (fig. 7.6, 7.7, and 7.8), 
exhibited less stability than expected. Also, system 2 was 
less stable than system 1. The largest difference in 
reading between the two receivers occurred at 0750 (fig. 
7.7). This difference may have been caused by the location 
of the antenna of system 2 under the antenna of system 1. 
At approximately 0800, the satellite constellation elevation 
angles were the largest, between 56 ° to 80° above the 
horizon (fig. 7.9), and antenna interference may have led to 
the large differences. Neither the satellite constellation 
azimuth angles (fig. 7.10) nor the user range errors (fig. 
7.11) appeared to have any effect on the stability. 

Comparison of the stationary GPS receiver positions with 
the position obtained with the geoceiver is depicted in 
figure 7. 12 for the stability and static tests. The 
comparison indicates that there is a nearly constant offset 
of the two positions with the GPS position almost due south 
of the geoceiver position. Further comparisons between the 
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stationary GPS receiver position and the geoceiver position 
reinforce this conclusion. In figures 7.13, 7.14, and 7.15, 
the mean distances between the two positions are 127, 173 
and 143 meters, respectively. The higher means for the data 
shown in figures 7.14 and 7.15 are due to the ORE (fig. 7.16 
and 7.17) for those dates and the NAVSTAR upload frequency. 
If the data from 13 May is deleted because of the high ORE 
for that date (fig. 7.16), and the median distance used for 
15 May is taken to be a representative value (130 meters on 
fig. 7.15), The difference is approximately 129 meters 
Between the two positions. The azimuth was also nearly 
constant with the observed median values of 181°, 207° and 
178 ° from south as shown in figures 7.18, 7.19, and 7.20 
respectively. Neglecting the data on 1 3 May because of ORE 
values (fig. 7.16), the remaining two values are nearly 

O 

identical and have a mean of 179.5 from south. Thus the 
data indicates a nearly constant offset between the position 
obtained with the geoceiver and the position obtained with 
GPS. This difference is about 130 meters and has an almost 
north-south orientation. 
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Figure 7.13 
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Figure 7.14 
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Figure 7.15 
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TOTAL USER RANGE ERROR 
13 HAY 1981 
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GREENWICH MEAN TIME (HOURS) 



Figure 7.17. 

TOTAL USER RANGE ERROR 
15 MAY 1981 
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GREENWICH MEAN TIME (HOURS) 



f RfcUUENC 1 £ S SAMPLE SUE 



Figure 7 . 18 . 
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F RE CUENC l E $ SAHPLE SIZE 



Figure 7.19 
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FPECUENCIES SAMPLE SIZE 



Figure 7.20 
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B. EVALUATION OF NONDIFFERENTIAL GPS 



Evaluation of the non-differential GPS positions was 
done on a set by set basis. Figures 7.21 through 7.31 depict 
the distance in meters between the theodolite-computed ship 
position and the shipboard GPS receiver position. Sets 1, 2 
and 3 have the smallest distance difference, as shown in 
figures 7.21, 7.22 , and 7.23, with mean values between 118 
and 133 meters and standard deviations ranging from 7.45 to 
16 meters. These smaller mean distances are a result of the 
low user range errors for 12 May as shown on figure 7.32. 
The highest mean distances occurred during sets 4, 9, 10, 
12, and 13 when the ORE values were high and the NAVSTAR 
uploads occurred (fig. 7.16 and 7.17). The remaining 
distances vary between 144 and 278 meters. 

The results during the best satellite observation 
periods are depicted in figures 7.33, 7.34, and 7.35, where 
the non-differential GPS ship position was compared to the 
theodolite-computed ship position, the distance between the 
two being plotted as a histogram. Mean distances on figures 
7.33 and 7.35 are representative of the five satellites 
observed and the improved accuracy of the system, with mean 
distances of 134 and 149 meters and standard deviations of 
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FRECUENC I £5 SAMPLE SIZE 



Figure 7.21. 
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FREQUENCIES SAMPLE SUE 



Figure 7>2 2 
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FREQUENCIES SAMPLE SUE 



Figure 7.23 
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FREQUENCIES SAMPLE SIZE 



Figure 7.24 
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F8ECUENCIES SAHPLE SIZE 



Figure 7.25 
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FREQUENCIES SAMPLE SIZE 



Figure 7.26 
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f*EQU*NClES SAMPLE SIZE 



Figure 7.27 
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FREQUENCIES SAMPLE SIZE 



Figure 7.28 
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FREQUENCIES SAMPLE SIZE 



Figure 7.29 
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FREQUENCIES SAMPLE SUE 



Figure 7.30 
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mcutNCics sample size 



Figure 7.31 
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Figure 7.32. 

TOTAL USER RANGE ERROR 
12 MAY 1981 
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GREENWICH MEAN TIME (HOURS) 



FPECUENCIES SAMPLE SIZE 



Figure 7.33 
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FRECUENC I E S SAMPLE SIZE 



Figure 7.34 
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F»cCJ£NC!£S SAHPLE SIZE 



Ffqure 7.35 
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10 and 52 meters. The large standard deviation for 15 May 
is apparent in the bimodal empirical frequency curve (fig. 
7.35) and is caused by the NAV STAS upload (fig. 7.17). in 
this case, the GPS receiver did not obtain the new ephemeris 
until halfway through the time period of the observations. 
The data for 13 May (fig. 7.34) reflects frequent NAVSTAR 1 
uploads and the resulting high ORB values (fig. 7.16) for 
that date. 

C. DIFFERENTIAL GPS EVALUATION WITH STATIC CORRECTION 

The results of the evaluation of the data of 12 May 
(sets 1, 2, and 3), computed using the differential GPS 
positions and the static corrections, are shown in figures 
7.36, 7.37, and 7.38. Mean difference between positions for 
the theodolite ship position and the GPS position are 29.1, 
27.4, and 64.8 meters, with standard deviations of 7.7, 7.5, 
and 22.7 meters. These distances are larger than those 
determined for the theodolite ship position and the 
differential GPS position (fig. 7.39, 7.40 and 7.41) . This 
higher mean distance indicates that the static correction, 
as applied, results in a larger error in the differential 
GPS position. Therefore, further evaluation of differential 
GPS with static corrections were not conducted. 
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f PECUENC I E S SAMPLE SIZE 



Figure 7.3 i 
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MecucNCiEs sample sue 



Figure 7.37 



» * • * • u.* • < 



O I 
I 
I 

♦ 

°: 

i 

♦ 

°i 

i 

♦ 
— i 
i 
i 

♦ 
o I 
i 

i 

psi I 
I 
I 

♦ 
*n i 
i 

i 

I 

I 

I 

+ 

O | 



rtu.tTttiittniuitmtmtmtirxi 7 :rrrryxrryrr rrrr 



o 

►— 

o 



z>- 



co^ooco 

♦ ♦♦♦♦ + * 

i i j UJ UJ 11 niiimii 

OOV'OC'C — 
j'rvj 

-r — r- '*■»©« 
y OT O 

©tnoao— m<*> 

— ._.< v»<NioO »*• 



— z— 

UJ <7UJ 

O — O 
ZQZ 

lit 



Z UJUJUJUJUJ 



X> ZZZZZ 
—0330000 

acrooooot 

• — « 

^ZO^'Otf'OX 

— — — •'‘ VUAf — c ^- < 

ox • • • • «x 



^ 

►-oooooo 

Z ♦> ♦ I I ♦ ♦ 

UJlU LU' u'jJujuj 
r ©cru>o— 
CNtr»^sr -ors» 
X *»f«l0©0* — 
OCyOcyao 

< fN or^PM'Oco 

3 t •••••• 

Z I III 
UJ 

O <✓»*✓> 

at jjvi 

UJ ZU-<N 

X 

O Ujac 

— Du/UJ 

T T T^*COCD 



O I 
I 



o i 
i 
i 



O I 
I 
I 

© I 
I 



^ X 

c— 

ex — 
2 
a. 

— z 
ia 

ui — 
►— 
«%J— • 

I u* 

►-a 

a. 

Z 

UJVl 

uj ex 

20 



'«-Z 

IUUJ 

*-nf 

UJ UJ 

TTu. 



— O 

u>< u 
IUI 

o — 

«?e 5 
i/l^t 
© 



—©— 0—0 

oooooo 

Uj UJUJUJUJUJ 

NNrg^ON 
C- C D O' — . 
oOcy^— © 
00 (T JtTiN 
Ovr^'NJtnaa 
N*T^<Drg^ 
• ••••• 
y \^ ryy % fr »«0 



o 

uj ac> *x 

o> -*uj uj 
QZUi>0 <X 
<T 'TO UJQ. 
uj— u.ZO'/* 
nc accuji ZO 
CX >r — Q| U <X — 

xx 



o in -r *— 
" j ‘'j'Tisr — 

~ vn n- j r o 
VO *' o^oon 

I"- U J. 



— — Z 



O ^ « <M <N nj ry o, ry 



UJ-t -1 
-J 7" ' a uj ‘ *j 

< r<>*zn 

*T 'tu'-'l 

— r t 

/•If 

Uj < • J V — — UJ *t 

->rr — r to r 



176 



fSECUINCIES SAMPLE SIZE 



Fioure 7.28 
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FREQUENCIES SAMPLE SIZE 



Figure 7. 39 
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FR 6 CUENCIES SAMPLE SIZE 



Figure 7.40 
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EAECUEftClES SAMPLE SIZE 
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D. EVALUATION OF DIFFERENTIAL GPS WITHOUT STATIC 
CORRECTIONS 

Evaluation of differential GPS without the static 
corrections were computed by sets for the five satellite 
constellations for each day of testing. Set 1 (fig. 7.39) 
shows a mean distance of 10.1 meters between the two 
positions and a standard deviation of 5.0 meters. The 
median azimuth between the two was 257° from south (fig. 
7.42) , which also was the approximate course of the vessel 
during the test. These values indicate that the error is 
caused by the bias, or alias, in GPS position due to the 
Manpack 24-second fixing interval and positioning algorithm. 
For this set, the GPS position was obtained after the 
theodolite position. 

During set 2 the mean distance difference between the 
two positions was 11.7 meters, and the standard deviation 
was 6.8 meters (fig. 7.40) ; the median azimuth was 71° 
(fig. 7.43) . These values were slightly greater than those 
of set 1 due to the increased URE for NAVSTAR 1 during the 
set (see 0700 on fig. 7.32) Additionally, the azimuth was 
nearly 180 ° from that computed in set 1, with the vessel on 
the same track as set 1. Unlike the error obtained from the 
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FPECUENCIES SAMPLE SUE 



Figure 7.42 
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FfteoutNCiES sample size 



Figure 7.43 
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Manpack bias on set 1 , the error in set 2 was due to a 
similar time delay bias because the theodolite position was 
obtained after the GPS position. 

Set 3 values on figures 7.41 and 7.44 are similar to 
those of set 2, except the mean distance values were 
greater — 44.7 meters for set 3 versus 11.7 meters for set 2. 
Mean distance increases for set 3 were due to the upload of 
NAVSTAR 1 (fig. 7.32) and the availability of only a four 
satellite constellation. Set 2 had a five satellite 
constellation available. 

For set 4, figures 7.45 and 7.46 show the distances and 
azimuths for the set. These histograms fail to illustrate 
the constantly increasing distances between the two 
positions. These increasing distances result from the 
uploads of NA7STAR' s 5 and 6 and the large ORE for NAVSTAR 1 
(fig. 7. 16) . Onboard the vessel the new ephemeris data 
uploaded for NAVSTAR *s 5 and 6 was not obtained until late 
in the set. Also near the end of the set, NAVSTAR 1 was 



uploaded with t 


he result 


that 


the last several 


distances 


were very large. 


as shown 


in f igure 7.16. 




In set 5, 


a very 


high 


ORE value for 


NAVSTAR 1 


contributed to 


the extremely 


large position 


difference 
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mCUENClES SAMPLE SIZE 



Figure 7.44 
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FLUENCIES sample size 



Figure 7.45 
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EKECUfNCIES SAMPLE SIZE 



Figure 7.46 
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between that determined by the theodolites and that for the 
differential GPS position (fig. 7.47). The azimuth remained 
nearly constant between the two positions with a median 
azimuth of 43° (fig. 7.48). These values for set 5 show 
that the large ORE for NA7STAR 1 was the most important 
contributing error source (fig. 7.16). The large ORE may 
have been caused by activation of the hydrazine 
stabilization jets to dampen momentum on NA7STAR 1. 

For set 6, the mean distance difference between the two 

O 

positions was 60.0 meters and a median azimuth of 93 (fig. 
7.49 and 7.50). These values indicate that the upload of 
NAVSTAR 1 near the beginning of the set (0715 on fig. 7.16) 
was received and utilized in the Hanpack position algorithm, 
but the higher distance between the positions also indicates 
a large ORE (also shown on fig. 7.16). 

Data from set 7 were not evaluated because the 
stationary Manpack receiver did not record the GPS position 
for each fix due to radio failure. The mean distance 
between the two positions for set 8 was 121 meters with a 
mean azimuth of 63 ° (fig. 7.51 and 7.52). Both the 
distances and azimuths show that the URE for NAVSTAR 1 at 
0530 (fig. 7.17) was increasing in magnitude and thus 
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FREQUENCIES SAMPLE SIZE 






Figure 7.47. 
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FREQUENCIES SAMPLE SUE 



Figure 7.48 



O 

I 

I 

♦ 

°l 

i 

O I 



o 

o 



o 

' 

o 



CO-— 

0000000 
♦♦♦♦♦♦♦ 
JJLUUJUjLU LUUJ 
— O — vriN|U^f«» 

OOirv— 0 >r 0 
ocooo'irv'Otn 
oocifimoo' 
r-^-J-lAON 
inofMrn<\i'TiA 
• •••••• 

—ii>rN»>rir\oo 



2 

a 



u. o 
u.— 
u. rsi 



O ZZ22Z 

3D r < <X <Z «* < 2 

iroooaot 

u^ZOinOiAOX 

— r»» o < 

ox 



(✓3 

>-000000 
2 ♦ ♦ ♦ I ♦ ♦ 

inininuJ n in 1 LU 

i r -<vO<oO 
Ococg— <\ iunO 
i/> ZO' s -r*0'r^'J > 

a. <j-n®u*« 
O J-OC'Nf-mO* 
►- cmO^'TO*-* 
o «•••••• 

UJ un®,4^INKD 
UC 2 I 



O 

o 



V»— 

UiV» 



— 2 
>- O iucr 
< — tDuJLU 

u> 



O 

- 

o 



U. — O 

a. O >rcuO-*cu-* 

u. X OOOOOO 



ztzitzitiizirzztztizzzirziiizzttrtttzzizizu.zzzzirr-o 

LU 1 -* 

u. 1 

LU 1 



z— 

O^C 

~Z 

— o 
(/>— 
CH 
a.— 



mu in nu n ml 

u>un — — r- wr\ 

i/>fMrg<jrarv* 
r-mfM<vi^>rg 
t*> or~ vn o 



«/i — — — 



CUO 
— Q. 
X 

\Z»V» 






o 

< 



NO 



UJ 0£> 

o>-*uj 
QZuj>q ac 
1 <r<xo ujq. 

H-—J lu— iLZO^I 
*t atfQiLi<2Q 
uj— a<v-cuj<i— • 
xh w»>«/ 3 oi:ac:x 



• -TU. 
UJO 



OOOOOOO 



LUUJ 1 1 in 11 1 11 1 hi 1 

33 -\j O' wn mco 

LU O— -OOvT 'f -O 

l/' X n O' O vn r>+~. 
X>- 3-^-i/w* 1 -*-^*-* 

»- uo^-rM^^ 
OO 

ZZ 

— «* Q— vr >r— «>r ^ 



*- zz 

LU 1 ** < 
_J ZZOLULU 

<1 Z-t<ZXZ 

ft <'UUi< 

►-z-rrtzz 

Z<Q-OUUa 

UJUJlUCt— — LU ■< 

uirntzoi 



190 



frecuejscies sample size 



Figure 7.49 
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FHECUENC1ES SAMPLE SUE 



Figure 7.50 
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32 1 S 3T«JKT$ ^31^3033^3 



Figure 7.51 
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FREQUENCIES SAMPLE SUE 



Figure 7.52 
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resulting in a greater positional difference. Another 
factor was the bias introduced by the Manpack 24-second 
reception period. This is evident in the azimuth for the 
theodolite ship position obtained from the GPS shipboard 
receiver. 

The data for sets 9 and 10 were combined for evaluation 
and the results were a mean distance between positions of 
241 meters and an azimuth of 267° (fig. 7.53 and 7.54). 
7alues of that magnitude show that the ORE for NA7STAR 1 was 
the error source. Data were sampled from 0600 to 0640 (fig. 
7.17). Thus the correlation between the high NA7STAR 1 ORE 
and the high mean distance difference between positions. 
Again, receiver bias also introduced an error into the 
evaluation. 

Analysis of set 11 shows a bimodal distribution in both 
the two position differences and azimuths (fig. 7.55 and 
7.56). The peaks on the distance histogram are 
approximately 133 and 160 meters; the peaks on the azimuth 
histogram, approximately 220 0 and 236° . Both peaks are 
caused by the upload of NA7STAR 3 halfway through the set 
and the increasing magnitude of the NAVSTAR 1 ORE just after 
the upload as shown by figure 7.17. The total positional 
difference for set 11 was approximately 154 meters. 
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F8ECUENC1ES SAMPLE SIZE 



Figure 7.53 
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f RECU6NC I ES SAMPLE SIZE 



Figure 7.54 
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FAECUENCIES SAMPLE SUE 



Figure 7.55 
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FREQUENCIES SAMPLE SIZE 



Figure 7.56 
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Evaluation of the combined sets 12 and 13 (fig- 7.57 and 
7.58) result in a mean distance between the two positions of 
172 meters and a mean azimuth of 262° . The mean distance 
was large because of the large magnitude of the NAVSTAR 1 
ORE (fig. 7. 17) . A contributing source of error was the 
bias introduced by the Manpack receiver. 

Overall, the data using a five-satellite constellation 
were as good or better than the data using only four 
satellites. The results for 12 May using five satellites 
(fig. 7.59) show a mean distance of 11.4 meters and a median 
of 10.6 meters. This compares favorably with the results 
from sets 1 and 2 (fig. 7.39 and 7.40), which had mean 
distances of 10.1 and 11.7 meters, respectively. The five 
satellite data from 13 May showed a very large variance due 
to the large NA7STAR 1 ORE. Therefore, data were not 
compared to remaining 13 May data. Data on 15 May for the 
five satellite constellation had a bimodal distribution 
because NAVSTAR 1 was uploaded. Once the uploading 
occurred, the peak distance was 137 meters (fig. 7.60) which 
compares favorably with 121, 241, 154, and 172 meters for 
sets 8, 9 and 10, 11, 12 and 13. Thus with a five-satellite 
constellation, positioning results were as good or better 
than the four-satellite constellation. 
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F8ECUENC 1 ES SAMPLE SIZE 



Figure 7.57 
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FREQUENCIES SAMPLE SIZE 



Figure 7.58 
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fPECUfNCIES SAMPLE SUE 



Figure 7.59 
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mUHNCieS SAHPLE SIZE 



Figure 7.60 
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Figures 7.61 to 7.74 show the difference by fix between 
the latitude and longitude plots of the theodolite ship 
position, the differential GPS position, and the uncorrected 
GPS position. During most of the testing, the longitude of 
both methods of GPS positioning was in closer agreement with 
the theodolite ship position than was the latitude. Shown 
on figure 7.9 and 7.10 are the satellite elevation and 
azimuth data for the testing period. This data indicates 
that the satellite constellation had a larger effect on 
latitude than longitude because the predominant ranging was 
computed from one quadrant. 

The large differences between the GPS position relative 
to the theodolite ship position were easily identified 
during the shipboard testing period. Latitude and longitude 
plots (fig. 7.65, 7.66, 7.69 through 7.74) for sets 4, 5, 9, 
10, and 12 show the deviations of the GPS receiver positions 
from the predicted positions that were obvious to the GPS 
operators during the shipboard testing. For sets 4 and 5, 
the stationary GPS receiver positions began to deviate 
greatly from the geoceiver position, resulting in a large 
differential correction to the GPS ship position. These 
large differential corrections caused the differential GPS 
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Figure 7.63. 
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Figure 7.65. 
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position shown on figures 7.65 and 7.66 to 



deviate a large 



amount from the theodolite position and the uncorrected GPS 
position. 
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VIII. CONCLO SIONS 



The objective of the study was to determine whether the 
application of GPS in a translocated or differential mode 
provides sufficient position accuracy for near-shore 
hydrographic surveys. Several assumptions were made prior 
to the test: (1) two GPS receivers required for the 
differential mode performed the computations similarly, so 
that in the worst case a constant bias would be displayed if 
both receivers were colocated; (2) the process of updating 
the ephemeris information would not contribute a significant 
positional error to GPS in the differential mode; (3) 
reasonable accuracy could be obtained by using the output 
from the receiver control display unit, that is, 
sophisticated processing methods were not required; and (4) 
the Hanpack receiver results would not differ when used in 
stationary or low-dynamic environments (shipboard at 4 
knots) . 

None of these assumptions proved to be valid. It was 
originally expected that the two receivers would produce 
similar positions if both sets were operating at the same 
time over the same location. One limitation of these 
receivers, however, is that current ephemeris information is 
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not received unless the receiver has been activated prior to 
the update. [Ref. 66] For a period of time after the 
receiver is activated, the portion of the satellite message 
that contains the updated ephemeris is not acknowledged. 
Therefore, position degrades rapidly when a NAVSTAR is 
updated and the current ephemeris information is not 
utilized by the receiver. For this study, unfortunately, 
the initial times of activation of the two receivers were 
not recorded. Hence, the effect on position caused by this 
receiver parameter was not defined. 

The ORE plots illustrate that the ephemeris updates did 
contribute significantly to the range error and, 
consequently, the resulting position error. The position 
determined by the Manpack receivers was not affected as 
seriously by the clock error in NAVSTAR 1 when a 
five-satellite constellation was used. The inaccuracy of 
pseudo- ranges received from this satellite required frequent 
updates from Tandenburg, especially on 13 Hay. A receiver 
that had the option of totally eliminating the erroneous 
pseudo-range would have given the best results for the 
study. Since the receivers operated independently, and the 
ephemeris information was not concurrent, the differential 
mode was effective for only a few cases. 
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A simple correction was applied to the values output on 
the control display unit of the Manpack receiver. This did 
not permit application of tropospheric or refraction error 
corrections, other than that applied by the receiver during 
computation of the position. A more sophisticated 
processing method, resulting in greater accuracy, would 
require acquisition of the pseudo- ranges and subsequent 
postprocessing to determine the receiver positions. 

The best results during the onboard portion of the study 
were obtained when the ship was allowed to drift at 
approximately 1 knot on 12 May. 

In conclusion, use of two Phase I Manpack receivers in 
the differential mode meets the required accuracy of 10 
meters only if the ship speed is less than 1 knot and other 
parameters such as ephemeris information and the number of 
satellites used for the position computation introduce a 
minimum amount of error. Further testing using differential 
GPS should incorporate more sophisticated processing methods 
and more advanced receivers to achieve the survey accuracy 
required for nearshore hydrographic surveys. 
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APPENDIX As. DESCRIPTION OF ERROR THEORY 



From application of error theory in the evaluation of 
positional information, it is possible to establish a 
meaningful accuracy statement subject to uniform 
interpretation. To provide a logical and acceptable basis 
FOB COMPUTATION AND COMPARISON, POSITIONAL ERRORS ARE 
assumed to follow a normal distribution. The assumption is 
valid because positional error components generally follow a 
normal distribution pattern when sufficient data is 
available. 

The statistical treatment of errors utilizes measurable 
errors obtained in the sources of positioning information. 
Analysis of the linear components provides a two-dimensional 
expression of the accuracy of the positioning system in the 
form of an error ellipse. The use of an error ellipse is 

i 

complicated by the problem of axes orientation and 
propagation of elliptical errors. Therefore, the ellipse is 
commonly replaced by circular fora which is easier to use 
and understand. The linear standard errors are combined and 
converted to a circular distribution and the final accuracy 
statement is expressed in terms of circular errors. 
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The circular form assumes the circular probability 



distribution function. This function expresses the 

probability that the radial error will be equal to or less 
than the radius of the probability circle. The dispersion 
of errors within the distribution is measured in terms of 
precision indexes, each of which represents the error which 
is unlikely to be exceeded for a given probability. The 

preferred circular precision indexes are the circular 
standard error (a ) , the circular probable error (CPE or 
CEP), the circular map accuracy standard (CMAS) , and the 
circular near- certainty error, or 3.5 sigma (3.5 o ). For 
each precision index, the given probability is 39.35%, 50%, 

90%, and 99.78%, respectively. 

Each of the precision indexes can be converted to one of 
the others by using the following table: 

Circular Error Conversion Factors 





39.35% 


50% 


90% 


99.78% 


39.35% 


1.0000 


1 . 1774 


2.1460 


3.5000 


50% 


0. 8493 


1.0000 


1.3227 


2.9726 


90% 


0.4660 


0.5486 


1.0000 


1.6309 


99.78% 


0.2857 


0.3364 


0.6131 


1 .0000 
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The rapid appro ximation for determining the circular 



standard error is 

CSE = 0.5000 (ax + ay) 

where a * and o y are linear standard errors obtained from the 
sources of positioning information, and therefore are the 
axes of the error ellipse. For a truly circular 

distribution, where °c is identical to the circular standard 
error, a x and a y are equal and the angle of intersection 
between the two is 90°. In all other cases, a normal 
circular error distribution is substituted for the 
elliptical distribution. The substitution is satisfactory 
for an analysis within specifiedcmin/omax ratios, where a min 
is the minor or smaller linear standard error of the two. 
Distortion in the error distribution between elliptical and 
circular forms at lower a min/ a max ratios limits the use of 
the circular concept to ratios greater than 0.2 . 

The CM AS is the precision index used as the O.S. 
National Map Accuracy Standard. It is interpreted as 
limiting the size of error which 90% of the positions will 
not exceed. CMAS is computed from circular standard error 
as follows: 

CMAS = 2 . 146 0 a c 
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The two-dimensional error distributions for positions 
determined by the ARTEMIS system and the theodolite network 
are described by CMAS. The accuracy of the specific 
positions is expressed by a statement of probability and 
error magnitude. The accuracy statement does not mean that 
the error in position is exactly the value shown in figure 
5.1, rather it expresses the probability that the error in 
position will not be larger than the error given. 
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